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CHAPTER  1 


INTRODUCTION 


The  loading  on  buried  structures  subjected  to  nuclear  or 


conventional  high  explosive  weapons  environments^  is  strongly  influenced 


by  the  backfill  adjacent  to  and  overlying  the  structure.  The  relative 


stiffness  of  the  structure  and  the  surrounding  soil  and  their 


interactions  will  determine  the  level  and  extent  of  damage  due  to  blast  sen 


Lbl<  Js 


loadings.  The  response  of  soil  to  high  amplitude,  short  duration, 


impulse  loads  is  thus  an  important  aspect  of  the  protective  construction 


design  and  analysis  (Calhoun  and  Kraft,  1966;  Schindler,  1968).  '■>  Soil 


response  may  be  investigated  in  the  laboratory^  -^However,  the  laboratory 


environment  must  be  able  to  reflect  the  type  of  confinement,  magnitude 


of  stress  change,  and  the  time  scale  of  loading  expected  in  the  problem  r 


(Whitman,  1970). 


A  reasonable  assumption  is  that  the  soil  response  in  the 


superseismic  region  of  an  atrblast  will  be  approximately  one-dimensional 


(Crawford,  Higgins,  and  Bultmann,  197*0.  The  time  scale  of  loading  is 


essentially  instantaneous  (Crawford  et  al.,  197*0.  Peak  stresses  can  be 


in  excess  of  1000  MPa  (Brode,  198*0.  To  model  the  soil  response  to  this 


type  of  loading,  uniaxial  strain  devices  have  been  developed  (e.g., 


Schindler,  1968).  Current  devices  can  apply  stresses  up  to  *100  MPa  with 


loading  times  on  the  order  of  0.3  milliseconds  (e.g.,  Jackson,  Ehrgott, 


and  Rohani ,  1 980) , 


The  split-Hopkinson  pressure  bar  (SHPB)  technique  has  been  used 


-  1  - 


m 
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over  the  last  thirty-five  years  as  a  tool  for  investigating  the  response 
of  metals,  rocks,  ceramics,  foams,  and  other  materials  to  short  duration 
compressive  impulse  loads  (e.g.,  Lindholm,  1964;  Hodge  and  Wasley,  1969; 
Christensen,  Swanson,  and  Brown,  1972).  Some  of  the  SHPB  devices  in  use 
can  apply  stresses  in  excess  of  1000  MPa  with  loading  times  on  the  order 
of  0.04  milliseconds  (Gaffney  and  Brown,  1984).  Until  recently,  the 
SHPB  technique  has  not  had  significant  application  in  the  field  of  soil 
mechanics.  Because  of  the  ability  to  apply  high  stresses  at  a  high  rate 
of  loading,  the  adaption  of  the  SHPB  technique  to  measure  the  dynamic 
response  of  soil  appears  to  be  a  promising  research  tool.  The  objective 
of  the  research  reported  herein  has  been  to  determine  whether  the  SHPB 
technique  can  be  effectively  used  to  measure  the  dynamic  response  of 
compacted  soil. 

The  use  of  soil  as  specimens  in  a  SHPB  experiment  is  not  a  trivial 
matter  because  soils  have  very  low  wave  speeds  (*  300  m/s)  in  comparison 
to  the  traditional  materials  tested  in  the  SHPB  (e.g.,  steel,  5000 
m/s).  Soils,  also  exhibit  nonlinear  hysteretic  behavior  which  will 
cause  the  amplitude  of  a  stress  wave  to  attenuate  as  it  propagates 
through  it  (Hendron  and  Auld,  1968).  In  addition,  the  relatively  low 
unconfined  compressive  strength  of  the  soil  (e.g.,  <  0.1  MPa)  creates 
difficulties  in  controlling  boundary  conditions. 

This  report  is  organized  into  seven  chapters.  Chapter  2  presents 
the  development  of  pressure  bar  experimentation  and  the  basic  theory  of 
measurement  used  in  reducing  the  experimental  data.  Chapter  3  describes 
the  SHPB  apparatus,  the  characteristics  of  the  soil,  and  how  the 
specimens  were  prepared.  The  data  reduction  procedure  have  been 


addressed  in  chapter  4.  Chapter  5  presents  the  experimental  results  in 


terms  of  the  soil  stress-strain  response,  and  establishes  that  the 
response  is  nearly  one  of  uniaxial  strain.  Experimental  replication  is 
also  demonstrated  and  the  observed  soil  behavior  at  different  compaction 
conditions  to  a  range  of  applied  stresses  are  presented.  Chapter  6 
discusses  the  results  presented  in  chapter  5  to  evaluate  the  observed 
stress-strain  response.  Strain-rate  independence  is  established  for 
specimen  strains  less  than  the  initial  gas  porosity.  In  addition, 
specimen  conditions  that  may  lead  to  erroneous  stress-strain  response 
and  an  evaluation  of  the  SHPB  experiment  are  presented.  Chapter  7 
presents  the  conclusions  drawn  from  the  analyses  and  provides 
recommendations  for  further  research.  The  appendices  supplement  and 
support  the  discussions  In  the  main  text. 


CHAPTER  2 


DEVELOPMENT  OF  THE  SPLIT-HOPKINSON 
PRESSURE  BAR  EXPERIMENT 

The  Hopklnson  Pressure  Bar 

In  191^,  Bertram  Hopklnson  devised  a  method  to  experimentally  study 
the  mechanical  behavior  of  solids  in  response  to  short  duration  impulse 
loads.  The  main  components  of  the  apparatus  constructed  by  Hopklnson 
were  a  long  cylindrical  steel  bar,  a  time  piece,  and  a  ballistic 
pendulum  (see  figure  2.1).  As  shown  in  figure  2.1,  the  cylindrical 
steel  bar  was  suspended  by  threads,  such  that  it  was  free  to  swing  in  a 
vertical  plane  and  remain  parallel  to  its  original  position.  The  time 
piece  was  the  same  diameter  as  the  cylindrical  bar  and  attached  to  it  by 
magnetic  attraction. 

The  experiment  was  initiated  by  the  impact  of  a  lead  bullet  or  the 
detonation  of  gun  cotton,  near  the  end  of  the  cylindrical  bar  opposite 
the  end  to  which  the  time  piece  was  attached.  The  resulting  impact 
produced  a  compressive  stress  wave  which  was  assumed  to  be  uniform  over 
the  cross  section  of  the  bar  and  to  propagate  down  the  length  of  the  bar 
and  through  the  joint  between  the  bar  and  time  piece  without 
distortion.  Upon  reaching  the  end  of  the  time  piece,  the  compressive 
stress  wave  was  assumed  to  reflect  as  a  tensile  wave.  Since  the  joint 
between  the  bar  and  the  time  piece  could  not  sustain  a  tensile  force 
greater  than  the  magnetic  force  connecting  them  the  time  piece  would 
separate  from  the  bar  when  the  tensile  force  across  the  joint  exceeded 
the  magnetic  force.  The  momentum  of  the  time  piece  was  measured  by  the 
ballistic  pendulum  and  assumed  to  be  equivalent  to  the  area  under  the 


stress-time  curve  between  points  of  equal  pressure,  separated  by  the 
time  required  for  the  wave  to  travel  the  length  of  the  time  piece 
twice.  By  measuring  the  variation  of  momentum  delivered  to  different 
length  time  pieces,  the  area  under  the  stress-time  curve  for  different 
time  intervals  could  be  obtained  allowing  a  complete  stress-time  curve 
to  be  constructed. 

Despite  its  innovative  approach,  Hopkinson's  method  suffered  from 
several  experimental  and  theoretical  limitations.  Experimentally,  the 
force  necessary  to  separate  the  bar  and  time  piece  was  unknown.  Also, 
the  available  instrumentation  could  not  detect  the  commencement  of  the 
time  interval  over  which  the  momentum  was  measured.  Theoretically,  it 
was  assumed  that  the  applied  stress  would  be  uniform  over  the  cross 
section  of  the  bar  and  that  the  stress  wave  would  propagate  down  the  bar 
without  distortion. 

These  limitations  were  addressed  and  overcome  through  an  extensive 
study  of  the  Hopkinson  bar  method  by  Davies  (19**8). 

Davies'  Contributions 

Davies  developed  a  means  of  measuring  the  displacements 
electrically  through  a  bar  condenser  unit  mounted  on  the  free  end  of  the 
pressure  bar.  Using  elastic  wave  theory,  Davies  translated  the 
electrical  records  into  a  precise  stress-time  curve.  The  electrical 
recording  system  avoided  the  use  of  a  time  piece,  and  alleviated  the 
associated  experimental  problems.  Therefore,  with  the  introduction  of 
electrical  data  recording,  most  of  the  experimental  limitations 
associated  with  Hopkinson's  original  experimental  method  are  avoided. 

As  part  of  the  same  study,  Davies  (19*18)  described  the  phenomenon 


of  wave  dispersion  and  established  the  accuracy  of  the  experimental 

results  when  assuming  one-dimensional  wave  propagation  in  the  pressure 

bar.  One- dimensional  elastic  wave  theory  assumes  that  a  wave 

propagating  along  a  bar  will  travel  without  change  of  form,  at  a 

constant  velocity  (C0),  that  is  independent  of  wavelength  (A).  The 

2 

constant  longitudinal  velocity  (CQ)  can  be  expressed  as  CQ  -  E/p,  where 
E  is  Young's  modulus  and  p  is  mass  density.  This  assumption  is  only 
true  if  the  wavelength  of  the  propagating  wave  is  long  compared  to  the 

radius  (R)  of  the  bar.  The  equation  governing  the  longitudinal 

vibrations  of  an  infinitely  long  circular  cylinder  (referred  to  in  the 
literature  as  the  dispersion  equation),  developed  independently  by 
Pochhammer  (1876)  and  Chree  (1889),  shows  that  the  one-dimensional 
assumption  will  only  be  approximate,  as  the  velocity  of  wave  propagation 
will  depend  on  A. 

A  compressive  stress  wave  generated  in  the  split  Hopkinson  pressure 
bar  experiment  is  composed  of  many  frequencies.  Because  the  wave 
velocity  (Cn)  depends  on  A,  each  frequency  component  will  travel  at  a 
different  Cn.  With  some  frequency  components  travelling  faster  than 
others,  the  wave  will  change  form  or  disperse  as  it  progagates.  The 

amount  a  wave  disperses  will  affect  the  uniformity  of  stress 

distribution  over  the  cross  section  of  the  bar.  The  more  a  wave 
disperses,  the  less  uniform  the  stress  distribution. 

By  using  the  dispersion  equation,  Davies  (19^8)  showed  that  the 
oscillations  in  the  displacement  history  recorded  at  the  surface  of  the 
bar  were  attributable  to  dispersion  and  not  experimental  inaccuracies. 
He  also  established  the  error  in  measuring  the  displacement  history  at 
the  surface  of  the  pressure  bar  to  be  ±  3  percent,  provided 


Kolsky* s  Contributions 

With  the  experimental  and  theoretical  foundations  of  the  method 
established,  Kolsky  (19^9)  modified  the  apparatus  to  permit  dynamic 
material  properties  to  be  determined  by  indirect  measurements.  Kolsky 
sandwiched  a  thin  cylindrical  specimen  (i.e.,  having  an  aspect  ratio, 
fc/d  -  0.1,  where  l  is  the  specimen  length,  and  d  is  the  specimen 
diameter)  between  two  pressure  bars.  The  pressure  bars  were  fitted  with 
condenser  microphones  for  data  recording.  He  then  applied  an  impact 
load  such  that  a  compressive  stress  wave  propagated  toward  the  specimen 
(see  figure  2.2).  When  the  propagating  stress  wave  reaches  the 

specimen,  a  portion  of  it  will  be  reflected  as  a  tensile  wave  and  a 
portion  will  be  transmitted  through  the  specimen.  The  amplitudes  of  the 
reflected  and  transmitted  waves  will  depend  on  the  physical  properties 
of  the  specimen  as  well  as  those  of  the  pressure  bars.  By  assuming  a 
uniform  distribution  of  stress  and  strain  along  the  longitudinal  axis  of 
the  specimen,  Kolsky  was  able  to  develop  relationships  between  the 
incident  and  transmitted  displacements  (recorded  at  the  respective 
condenser  microphones  mounted  on  the  bars)  and  the  average  stress, 
strain,  and  strain-rate  in  the  specimen.  The  derivation  of  these 
equations  is  presented  in  Appendix  B.  This  experimental  technique  is 
now  known  as  the  Kolsky  method,  or  the  split-Hopkinson  pressure  bar 
(SHPB)  method. 

By  using  thin  specimens,  Kolsky  attempted  to  avoid  the  complication 
of  axial  inertia.  Radial  inertia  effects  were  accounted  for  through  a 
numerical  correction  that  assumed  frictionless  specimen-bar 


interfaces 


The  effects  of  axial  and  radial  inertia  in  a  SHPB 


experiment  are  important  because  they  act  to  oppose  the  equilibration  of 
stress  within  a  specimen.  Kolsky  determined  that  by  using  a  thin  layer 
of  lubrication  between  the  specimen,  and  the  bars,  the  specimen-bar 
interfaces  could  be  assumed  to  be  frictionless. 

Contributions  of  Recent  Investigators 

Since  the  modification  made  by  Kolsky  to  the  original  Hopkinson  bar 
experimental  configuration,  the  only  significant  change  has  been  the 
introduction  of  strain  gauges  for  date  recording.  The  use  of  strain 
gauges  was  first  reported  by  Hauser,  Simmons,  and  Dorn  (1961). 
Replacement  of  the  condenser  microphones  with  strain  gauges  avoids  the 
requirement  of  differentiating  the  displacement  records  as  a  continuous 
strain  history  can  now  be  recorded. 

To  simplify  the  equations  for  stress  and  strain  (see  Appendix  B) 
some  investigators  have  assumed  the  stress  along  the  axis  of  the 
specimen  to  be  uniform  (e.g.,  Lindholm,  196*0.  This  assumption  was  not 
used  in  this  investigation.  Instead,  all  computations  of  specimen 
stress,  strain,  and  strain-rate  have  been  performed  using  equations 
developed  in  Appendix  B.  A  discussion  of  the  experimental  assumptions 
is  presented  in  appendix  D. 

A  critical  analysis  was  made  of  the  Kolsky  technique,  and  its 
assumptions  by  Davies  and  Hunter  (1963).  While  their  experimental 
technique  was  essentially  the  same  as  that  used  by  Kolsky,  by  using  an 
analytic  method,  a  criterion  was  developed  to  estimate  whether  the 
assumptions  of  the  experiment  are  satisfied,  principally  the  uniform 
stress  condition  (see  appendix  D). 

Davies  and  Hunter  (1963)  determined  that  to  minimize  the  effects  of 
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friction,  the  specimen  aspect  ratio  should  be  approximately  unity.  The 

criterion  used  to  estimate  when  stress  equilibrium  in  the  specimen  has 

been  achieved  is  based  on  the  duration  of  the  applied  stress  wave  (T) 

and  the  characteristic  velocity  at  which  the  wave  propagates  through  the 

2 

specimen.  The  characteristic  velocity  (C  -  E  /p  ,  where  E_  is  the 

3  S  S  ® 

slope  of  the  tangent  to  the  stress-strain  curve  at  a  point 

(da/de)  and  p  is  the  density  of  the  specimen)  is  that  given  by  the 
s 

Taylor-von  Karmen  theory  for  a  plastically  deforming  specimen  (von 
Karmen,  19*12;  Taylor,  1958).  Combining  energy  principles  with  the 
Taylor-von  Karmen  theory,  the  derived  criterion  for  estimating  when 
stress  equilibrium  in  the  specimen  has  been  achieved  is  given  as  (Davies 
and  Hunter,  1963); 


(2.1) 


Davies  and  Hunter  (1963)  suggested  that  when  this  inequality  is 
violated,  equilibrium  has  not  been  reached  and  the  stress-strain  data 
may  be  in  error. 

The  SHPB  experimental  environmental  was  also  studied  numerically  by 
Bertholf  and  Karnes  (197*0.  They  used  a  two-dimensional,  finite- 
difference,  elastic-plastic,  wave  propagation  computer  code  TOODY 
(Bertholf  &  Benzley,  1968)  to  investigate  the  response  of  an  aluminum 
specimen  which  was  modeled  as  a  strain-rate  independent  material.  They 
determined  that  the  SHPB  experiment  could  be  used  to  accurately,  and 
reliably  determine  material  response  at  high  rates  of  loading,  provided 
that  certain  precautions  were  taken  to  minimize  the  effects  of  inertia 
and  friction.  It  was  shown  that  the  effects  of  inertia  and  friction 
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between  the  specimen  and  the  bars  could  lead  to  the  conclusion  that  the 


» 


material  response  was  rate  dependent  (i.e.,  inertia  and  friction 
affected  the  degree  to  which  stress  and  strain  uniformity  in  a  specimen 
could  be  achieved),  especially  for  thf  case  of  relatively  thin  specimens 
(i.e.,  small  aspect  ratios).  Their  numerical  computations  showed  that 
the  influence  of  inertia  on  experimental  data  could  be  minimized  by 
bounding  the  maximum  specimen  strain-rate  such  that  (Bertholf  and 
Karnes,  1974); 

De  -  5  x  103  cm  s_1  ,  (2.2) 

max 

where  the  loading  wave  is  a  ramp  and; 

T 

>  16  ps  cm  1  ,  (2.3) 

where  tp  is  the  rise-time  of  the  loading  wave  and  D  is  the  diameter  of 
the  pressure  bars.  In  addition,  they  showed  that  if  sufficient  care  is 
taken  to  lubricate  the  ends  of  the  specimen-bar  interfaces,  the 
influence  of  friction  on  experimental  results  can  be  minimized.  Their 
investigation  also  determined  that  the  criterion  proposed  by  Davies  and 
Hunter  to  determine  if  the  reduced  data  may  be  in  error,  was  indeed 
reasonable. 

Since  Davies'  pioneering  work  several  investigators  have  throughly 
examined  the  dispersive  nature  of  wave  propagation  in  elastic  bars 
(e.g.,  Curtis,  I960;  Yeung  Wye  Kong,  Parsons,  and  Cole,  1974).  However, 
only  recently  has  a  numerical  procedure  to  account  for  wave  dispersion 
been  developed  that  can  easily  be  incorporated  into  the  standard  SHPB 
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data  reduction  technique  (Follansbee  and  Frantz,  1983;  Follansbee  and 
Frantz,  1984).  This  technique  has  been  adapted  for  use  in  this 
research. 

As  previously  mentioned  the  generated  stress  wave  is  not  composed 
of  a  single  frequency,  but  rather  a  spectrum  of  frequencies,  with  each 
frequency  travelling  at  its  own  respective  phase  velocity  (Cn)  and 
wavelength  (A).  Due  to  this  frequency  variation,  oscillations  develop 
in  the  propagating  stress  wave  which  have  no  counterpart  in  the  applied 
stress  wave  (Curtis,  I960).  The  nature  of  the  wave  oscillations  can  be 
described  mathematically  through  the  use  of  the  dispersion  equation. 
The  dispersion  equation  can  be  written  in  the  form  (Wasley,  1973); 

(h2  ♦  k2)  J , ( gR )  Jj(hR)  -  (h2  -  k2)*  J0(gR)  J,(hR)  (2.4) 
-  4  k2ghJ0  (hR)  J,(gR)  -  0  , 


where 
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Jq  is  a  Bessel  function  of  the  first  kind  of  order  zero,  J1  is  a  Bessel 
function  of  the  first  kind  of  order  one,  R  is  the  radius  of  the  bar,  k 
is  the  wave  number  (2ir/A),  A  and  y  are  Lame's  constants,  p  is  the 
circular  frequency,  and  p  is  the  bar  density.  The  solution  to  the  first 
mode  of  vibration  of  the  dispersion  equation  is  shown  in  figure  2.3  for 
a  material  with  a  Poisson’s  ratio  of  0.29. 

The  plot  of  the  fundamental  mode  of  the  dispersion  equation  shows 


that  a  high  frequency  wave  will  travel  slower  than  a  low  frequency  wave 
(see  figure  2.3).  Hence,  as  a  wave  propagates,  the  higher  frequency 
components  of  the  wave  will  lag  behind  the  lower  frequency  components 
causing  a  change  in  the  original  shape  of  the  wave.  This  wave 
dispersion  as  related  to  the  SHPB  experiment  increases  the  difficulty  in 
interpreting  specimen  behavior  as  deduced  from  the  standard  data 
reduction  technique  (Follansbee  and  Frantz,  1983). 

An  assumption  of  the  dispersion  correction  procedure  is  that  only 
the  fundamental  mode  of  vibration  is  excited  during  the  SHPB 
experiment.  This  assumption  is  supported  by  the  work  of  Davies  ( 19*48) 
and  Follansbee  and  Frantz  (1983).  Based  on  this  assumption,  the 
fundamental  mode  of  vibration  of  the  dispersion  equation  has  been 
approximated  by  fitting  the  data  with  an  equation  of  the  form 
(Follansbee  and  Frantz,  1983); 


jr-  =  0.576*4  + 
o 
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(2.5) 

Figure  2.  *4  compares  the  approximation  to  the  solution  of  the  dispersion 
equation  for  the  fundamental  mode  of  vibration.  Also,  the  applied 
stress  wave  has  been  represented  as  a  Fourier  cosine  series; 


f(t)  -  g-  +  l  Dn  cos(nu)Qt  -  $)  , 
n»1 


(2.6) 


where  AQ  is  the  amplitude  of  the  largest  frequency  component,  u>Q  is  the 
lowest  frequency  component,  Dn  is  the  amplitude  of  frequency 
component  nu)Q,  and  $  is  the  phase  angle.  The  phase  angle  is  the 


component  used  to  correct  the  wave  form  for  dispersion.  A  complete 
mathematical  description  of  the  dispersion  correction  procedure  is 
presented  in  appendix  A. 

Application  of  the  SHPB  to  Soils 

In  1967  Fletcher  and  Poorooshab  presented  a  paper  on  their  work 
with  the  SHPB  using  thin  clay  specimens.  The  purpose  of  their  study  was 
to  investigate  the  response  of  a  kaolin  clay  to  loads  of  low  magnitude, 
applied  at  high  rates.  The  largest  magnitude  load  was  less  than  1  MPa 
and  the  average  strain-rate  was  200  s.  The  loading  stress  wave  was 
initiated  by  a  steel  ball  accelerated  to  the  desired  impace  velocity  by 
rolling  it  down  a  ramp. 

Their  results  showed  that  the  response  of  the  clay  was  influenced 
by  the  preconsolidation  pressure  at  which  it  was  prepared.  At  applied 
loads  greater  that  the  preconsolidation  pressure,  a  peak  stress  was 
achieved  that  was  greater  than  the  preconsolidation  pressure.  After 
reaching  the  peak,  the  specimen  stress  fell  rapidly  to  a  level  at  which 
the  clay  appeared  to  flow  at  a  constant  stress.  For  applied  loads  less 
than  the  preconsolidation  pressure,  the  peak  stress  achieved  was  equal 
to  the  applied  stress.  Once  this  stress  level  was  reached,  the  clay 
flowed  at  that  stress  with  no  loss  of  strength,  in  contrast  to  the 
specimens  where  the  applied  stress  was  greater  than  the  preconsolidation 
pressure.  The  equations  they  present  to  determine  the  specimen  stress 
and  strain  assume  that  the  stress  is  uniform  across  the  specimen.  This 
assumption  is  questionable  for  soil  specimens. 

A  more  recent  investigation  into  the  response  of  soils  using  the 
SHPB  has  been  conducted  by  Gaffney,  Brown,  and  Felice  (1985).  The  soil 


was  a  clayey  sand  obtained  from  the  CARES-Dry  site,  located  on  Luke  Air 
Force  Base,  Arizona.  The  specimens  were  prepared  by  static  compaction 
in  thick-walled  confining  cylinders.  Four  specimens  were  prepared  at  a 
moisture  content  of  3.5  percent  and  at  a  wet  density  of  1.77  g/cm^.  Two 
nominal  specimen  lengths  were  used,  13  mm,  and  25  mm.  The  initial 
volume  of  air  contained  in  the  specimen  was  approximately  29  percent  of 
the  total  specimen  volume.  The  applied  stresses  for  the  experiments 
ranged  up  to  500  MPa. 

In  addition  to  the  SHPB  experiments,  gas  gun  experiments  at  high 
strain-rates  (>  5000  s_1 )  and  quasi-static  experiments  (strain-rates  on 
the  order  of  (5  X  10~3  s-1)  were  also  conducted.  Through  comparison  of 
the  results  obtained  over  a  range  of  strain-rate  regimes,  it  was 
concluded  that  for  strain-rates  below  5000  s'1  the  response  of  the  soil 
with  a  moisture  content  of  3-5  percent  was  independent  of  strain-rate. 
It  was  noted  that  in  the  SHPB  experiments,  the  specimens  retained  some 
volume  of  air  and  that  the  soil  response  may  be  significantly  different 
if  all  the  air  voids  had  been  closed  (Gaffney  et  al.,  1985).  The  same 
SHPB  apparatus  used  by  Gaffney  et  al.  (1985)  was  also  used  for  this 


investigation. 


figure  2.3  Solution  to  the  fundamental  mode  of  vibration  of  the  dispersion 
equation  (adapted  from  Davies,  1946) 


CHAPTER  3 


SOIL  SPECIMEN  PREPARATION  AND  EXPERIMENTAL  APPARATUS 

Soil  Characteristics 

The  soil  used  for  the  experiments  described  herein  was  sampled  in 
bulk  quantities  from  a  test  pit  in  an  undisturbed  area  of  the  McCormick 
Ranch  test  site  located  on  Kirtland  Air  Force  Base,  New  Mexico,  The 
site  has  been  used  for  explosive  testing  by  the  Department  of  Defense 
since  the  early  1960's.  The  geology  of  the  site  consists  of  some 
surface  playa  underlain  by  unconsolidated  eolian  sand,  alluvium,  and 
lacustrine  deposits  15  to  30  m  thick  (Bedsun,  1983).  The  water  table  is 
estimated  to  lie  at  a  depth  of  60  to  90  m  (Bedsun,  1983). 

In  order  that  the  soil  be  as  free  of  organic  material  as  possible, 
the  surface  vegetation  was  scraped  away  and  the  sample  taken  from  a 
depth  of  1  to  2  m.  After  arriving  at  the  University  of  Utah  soils 
laboratory,  the  samples  were  mixed  to  break  apart  large  clumps  and  to 
achieve  a  uniform  mixture.  The  following  tests  were  then  performed  to 
determine  fundamental  soil  properties:  (1)  mechanical  sieve  and 
hydrometer  analyses  to  determine  grain-size  distribution,  (2)  Atterberg 
limits,  (3)  hygroscopic  moisture  determination,  and  (4J)  specific 
gravity.  These  tests  were  conducted  in  accordance  with  the  American 
Society  for  Testing  and  Materials  standards  ( ASTM,  1982).  The 
moisture/density  relationship  was  determined  by  the  Harvard  miniature 
(kneading)  compaction  procedure  using  three  lifts  with  twenty-five  *40 
pound  t3mps  per  lift.  An  x-ray  diffraction  study  was  also  performed  to 
determine  the  soil  constitutive  minerals. 


The  results  of  these  tests  showed  the  soil  to  be  a  nonplastic, 


clayey  silty  sand  (see  figure  3.1)  with  a  specific  gravity  of  2.67  and 
hygroscopic  moisture  content  of  2.9  percent.  The  constitutive  minerals 
and  their  percentages  are,  quartz,  60  percent,  feldspar,  30  percent, 
with  traces  of  hornblende,  illite,  montmorillonite,  cristobalite,  and 
calcite.  A  carbonate  test  showed  the  percentage  of  calcite  to  be 
approximately  2  percent. 

Specimen  Preparation 

In  order  to  control  the  density  and  moisture  of  test  specimens, 
they  were  prepared  using  static  compaction  procedures  to  compact  a  known 
weight  of  soil  into  a  specified  volume.  The  majority  of  specimens  were 
prepared  near  the  optimum  moisture  content  (13.3%)  and  dry  density  (1.87 
g/cm^)  as  determined  by  the  Harvard  miniature  compaction  procedure  (see 
figure  3-2).  A  complete  tabulation  of  the  initial  specimen  parameters 
for  each  specimen  is  presented  in  appendix  E.  To  minimize  variability 
among  specimens,  the  soil  was  mixed  in  batches  sufficient  to  provide  all 
specimens  for  a  particular  test  session  (a  minimum  of  five  experimental 
specimens).  Before  the  addition  of  water,  the  soil  was  passed  through  a 
No.  4  sieve  (4.75  mm  opening),  weighed  and  placed  in  a  large  flat  pan. 
The  specified  amount  of  moisture  was  added  by  using  a  spray  bottle  so 
that  an  even  distribution  could  be  obtained.  The  soil  was  then  mixed 
throughly.  After  mixing,  a  damp  cloth  was  placed  over  the  soil  for  a 
period  of  20  minutes  to  allow  equilibration  of  the  soil-water  mixture. 
Following  the  equilibration  period,  the  soil  was  again  mixed  for 


uniformity  prior  to  compaction.  At  this  time  a  small  specimen  of  the 
batch  mix  was  taken  for  a  moisture  content  determination  using  the 
standard  laboratory  technique  for  moisture  content  determination  ( ASTM- 


Specimens  were  statically  compacted  in  the  confining  cylinders  by 
using  a  hydraulic  press  with  spacer  rings  to  control  the  specimen  length 
(i.e.,  density).  Two  nominal  specimen  lengths  (12.7  mm  and  6.35  mm) 
were  chosen  based  on  preliminary  tests  which  showed  that  for  specimen 
lengths  greater  than  13  mm  difficulty  was  encountered  in  obtaining  a 
transmitted  stress  signal.  Each  specimen  was  then  sealed  in  a  plastic 
bag  to  minimize  any  moisture  loss  that  might  occur  prior  to  the 
experiment.  The  same  preparation  procedure  was  followed  for  specimens 
prepared  on  the  wet  and  dry  side  of  optimum  conditions. 

Eight  specimens  were  prepared  for  the  sole  purpose  of  quantifying 
the  amount  of  moisture  loss  to  be  expected  between  the  time  the  specimen 
was  prepared  and  the  time  the  experiment  was  performed.  These  specimens 
were  prepared  in  a  equivalent  fashion  to  the  test  specimens.  However, 
after  placement  in  the  SHPB  apparatus,  these  specimens  were  removed  and 
the  moisture  content  of  the  entire  specimen  was  determined.  The  results 
of  these  experiments  showed  an  average  loss  of  5.4  percent  of  the  total 
moisture  with  a  standard  deviation  of  1.8  percent.  These  results  are 
believed  to  be  representative  of  the  pre-experiment  moisture  losses  for 
the  specimens  prepared  near  optimum  compaction  conditions.  The  results 
for  each  of  these  tests  may  be  found  in  appendix  E. 

Prior  to  performing  each  experiment,  the  specimen  was  removed  from 
the  plastic  bag  and  its  length  was  measured  inside  the  confining 
cylinder  by  lowering  a  depth  gauge  from  the  lip  of  the  cylinder  to  the 
surface  of  the  specimen  on  both  ends  of  the  cylinder.  At  the  completion 
of  each  experiment,  the  specimen  length  was  again  measured.  The 
specimen  was  then  extruded  from  the  confining  cylinder  and  a  portion  of 


the  specimen  used  to  determine  the  post  experiment  moisture  content. 

Whenever  a  specimen  of  soil  is  compacted  using  static  compaction 
procedures  there  will  exist  moisture  and  density  variations  over  its 
length.  A  limited  experimental  study  was  conducted  to  evaluate  the 
extent  of  these  variations  for  specimens  compacted  using  the  static 
compaction  technique  followed  in  this  study.  The  nominal  length  of  the 
specimens  prepared  to  study  these  variations  was  12  cm.  The 
experimental  procedure  and  results  are  presented  in  appendix  F.  While 
there  is  some  scatter  in  the  data,  the  results  indicate  that  in  general, 
a  specimen  will  have  the  highest  density  at  its  ends  and  the  lowest 
density  at  its  center.  The  moisture  content  will  vary  from  a  high  at 
the  center  to  lower  at  the  ends. 

Although  the  specimen  lengths  used  in  the  moisture  and  density 
variation  study  greatly  exceeded  the  specimen  lengths  used  for  the  SHPB 
experiments,  the  results  can  readily  be  used  to  estimate  the 
moisture/ density  variations  in  the  SHPB  specimens.  In  short  sections 
near  the  center  of  the  specimen,  the  variations  in  density  and  moisture 
were  less  than  2  and  4  percent;  respectively.  During  compaction,  these 
variations  are  affected  by  friction  between  the  soil  and  the  confining 
cylinder.  For  long  specimens  a  greater  area  will  exist  over  which 
friction  forces  can  act.  Hence,  friction  will  influence 

moisture/ density  variations  to  a  greater  degree  in  longer  specimens  than 
it  will  in  short.  Therefore,  these  values  can  be  considered  upper 
bounds  for  the  short  specimens  used  in  the  SHPB  experiments. 

Specimen  Confining  System 


To  achieve  a  nearly  uniaxial  strain  environment  for  the 


experiments,  the  soil  specimens  were  compacted  in  a  thick-walled 
confining  cylinder  which  helped  enforce  the  condition  of  uniaxial 
strain.  The  concept  of  confining  a  SHPB  specimen  such  that  it  will 
experience  a  condition  of  nearly  uniaxial  strain  was  first  presented  by 
Bhushan  and  Jashman  (1978). 

The  nominal  dimensions  of  the  confining  cylinder  are  60.3  mm  inside 
diameter,  102  mm  outside  diameter,  and  H4.5  mm  in  length.  The  confining 
cylinder  served  several  purposes;  first,  to  contain  the  soil  specimen 
itself,  and  second,  because  the  specimen  will  experience  a  state  of 
nearly  uniaxial  strain,  the  effects  of  radial  inertia  are  avoided.  The 
confinement  also  tended  to  prevent  specimen  distortion  or  barrelling 
during  the  experiment. 

When  the  confining  cylinder  containing  the  specimen  was  placed 
between  the  pressure  bars  approximately  19.0  mm  of  the  cylinder 
overlapped  the  bars  on  each  end  (see  figure  3-3).  To  determine  if  the 
confining  cylinder  was  transferring  any  stress  to  the  transmitter  bar,  a 
test  was  conducted  with  the  bars  separated  a  distance  of  3.0  mm  and  the 
confining  cylinder  placed  over  the  air  gap.  The  results  indicated  that 
the  confining  cylinder  did  not  transfer  any  measurable  stress  to  the 
transmitter  bar. 

Experimental  Apparatus 

The  Split-Hopkinson  pressure  bar  apparatus  used  in  this  research  is 
the  property  of  the  Geophysics  Group  at  the  Los  Alamos  National 
Laboratory,  Los  Alamos,  New  Mexico.  A  photograph  of  the  SHPB  apparatus 
is  shown  in  figure  3-  ^ »  while  figure  3.5  presents  a  diagrammatic 
representation  of  the  same  equipment.  The  main  components  of  the  system 
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are  the  gas  gun,  the  reaction  frame,  and  the  incident  and  transmitter 
pressure  bars  (additional  components  are  noted  in  figure  3.5).  Figure 
3-6  presents  an  enlarged  schematic  of  the  incident  and  transmitter 
pressure  bars  showing  the  specimen  positioning  and  the  location  of  the 
strain  gauges. 

The  incident  and  transmitter  pressure  bars  are  constructed  of 
Vascomax  350  CVM  maraging  steel  that  has  been  heat  treated  to  sustain  a 
yield  stress  of  about  2  GPa.  Each  pressure  bar  is  60.3  mm  in  diameter 
and  1.22  m  in  length.  The  bars  ride  on  adjustable  teflon  bearings  that 
allow  unrestricted  motion  in  the  horizontal  plane.  The  teflon  bearings 
do  not  restrict  the  passage  of  the  stress  wave  (Lindholm  1 964 ) . 
Adjustments  of  the  teflon  bearings  for  system  alignment  are  made  through 
four-wedge  supported  mounting  rings.  The  entire  system  is  supported  by 
a  four-column  reaction  frame  mounted  on  a  steel  I-beam. 

The  applied  stress  wave  is  initiated  by  the  impact  of  the  striker 
bar  on  the  incident  bar.  The  striker  bar  constructed  of  the  same 
material  has  a  slightly  larger  diameter  (60.5  mm)  than  the  pressure 
bars.  Three  striker  bar  lengths  were  available,  0.127  m,  0.251*  m,  and 
0.508  m.  To  minimize  misalignment  during  impact  and  increase  the  rise 
time  of  the  applied  stress  wave,  the  impact  end  of  the  striker  bar  is 
slightly  rounded  (Follansbee  and  Frantz,  1983).  The  striker  bar  is 
accelerated  down  the  launch  tube  to  the  appropriate  impact  velocity  by 
the  gas  gun.  The  amplitude  of  the  applied  stress  wave  is  proportional 
to  the  impact  velocity  of  the  striker  bar  which  is  controlled  by  the  gas 
gun  breech  pressure  (Nagy  and  Muelenhaupt,  1983).  Appendix  C  presents  a 
derivation  of  the  relationship  between  the  impact  velocity  of  the 
striker  bar  and  the  amplitude  of  the  applied  stress  wave.  The  end  of 


the  launch  tube  la  vented  so  that  the  gas  driving  the  striker  bar  is 
expelled  before  impact,  to  avoid  multiple  impacts  by  the  striker  bar  on 
the  incident  bar.  Just  prior  to  impact,  the  striker  bar  velocity  is 
measured  by  three  sets  of  diode  lasers  and  photodectectors  mounted  in 
the  end  of  the  launch  tube. 

The  applied  stress  wave  in  the  bars  is  monitored  by  resistance 
strain  gauges  mounted  on  the  radial  surface  of  the  pressure  bars  (see 
figure  3-6)  The  strain  gauges  are  mounted  in  pairs  on  opposite  sides  of 
the  respective  pressure  bars  and  connected  in  a  half-bridge 
configuration  to  nullify  any  bending  strains.  The  strain  gauges  are 
bonded  to  the  pressure  bars  at  mid'length.  The  positioning  of  the 
strain  gauges  on  the  incident  bar  is  particularly  critical  because  it  is 
necessary  to  obtain  continuous  records  of  the  incident  and  reflected 
waves  without  overlap. 

The  data  recorded  from  strain  gauge  bridges  are  filtered  and 
preamplified  and  then  routed  to  a  data  acquisition  system.  The  data 
were  then  transferred  by  magnetic  tape  to  a  CRAY-1S  mainframe  computer 
located  at  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 
New  Mexico  for  complete  data  processing. 
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CHAPTER  4 


DATA  ANALYSIS 


Experimental  Assumptions 

The  SHPB  experimental  method  is  limited  by  how  well  the  assumptions 
required  to  reduce  the  data  are  satisfied.  These  assumptions  are: 

(1)  there  is  a  uniform  distribution  of  axial  and  radial  stress 
over  the  length  of  the  specimen, 

(2)  the  waves  in  the  pressure  bars  propagate  without  dispersion, 

(3)  the  stress  state  over  the  cross  sectional  area  of  the  pressure 
bars  is  one-dimensional,  and 

(4)  the  interfaces  between  the  pressure  bars  and  the  specimen  are 
frictionless. 

These  are  the  basic  assumptions  made  Kolsky  (1949).  Each  assumption  and 
its  bearing  on  the  experimental  results  are  described  in  appendix  D. 

Data  Reduction  Procedure 

This  section  will  present  the  steps  used  in  the  data  reduction 
procedure  to  compute  the  average  specimen  strain-rate-time,  strain-time, 
and  stress-time  response. 

The  steps  are: 

(1)  data  input, 

(2)  correct  each  wave  for  dispersion, 

(3)  compute  the  average  specimen  strain-rate-time,  strain-time, 
and  stress-time  response  using  equations,  and 

(4)  output  processed  data. 

To  perform  the  computations  required  in  the  data  reduction 
procedure  a  computer  program  (REDUCE)  was  developed.  A  flowchart  and 
listing  of  REDUCE  is  presented  in  appendix  G  along  with  an  example  of 
the  output  for  experiment  134. 


Examples  of  the  raw  data  recorded  by  the  strain  gauges  for  the  two 
nominal  specimen  lengths  are  shown  in  figures  4.1  and  4.2.  Figures  4.3 
and  4.4  compare  the  raw  data  of  figures  4.1  and  4.2  with  the  data  after 
the  dispersion  correction  procedure  has  been  performed.  Because  the  raw 
data  are  recorded  at  some  distance  away  from  the  specimen-bar 
interfaces,  and  the  equations  used  in  the  computations  are  only  valid  at 
the  specimen-bar  interfaces,  each  wave  must  be  corrected  for  the 
dispersion  that  occurs  in  the  incident  and  transmtter  bars  between  the 
strain  gauge  and  the  specimen-bar  interface. 

Figures  4.5  and  4.6  illustrate  the  waves  of  figures  4.1  and  4.2 
after  correction  of  dispersion  and  time  shifting  to  the  specimen-bar 
interfaces.  The  incident  and  reflected  waves  are  shifted  to  the 
specimen- incident  bar  interface  so  that  they  will  start  at  the  same 
moment.  The  transmitted  wave  will  be  shifted  to  the  specimen- 

transmitter  bar  interface  but  with  a  delay  which  depends  on  the  wave 
velocity  in  the  specimen.  As  the  distance  between  the  specimen- bar 
interface  and  the  strain  gauge  is  the  same  for  both  pressure 
bars,  Ax,  each  wave  is  time  shifted  by  Ax/Co>  where  CQ  (4900  m/s)  is  the 
wave  velocity  in  the  pressure  bars. 

Examples  of  the  average  strain-rate-time,  strain-time,  stress-time, 
and  stress-strain  response  as  computed  by  REDUCE  for  the  two  nominal 
specimen  lengths  are  shown  in  figures  4.7  through  4.14.  Appendix  I 
presents  reduced  stress-strain  curves  for  all  specimens  tested. 


Uncorrected  and  corrected  stress  waves  for  a  6.35mm  specimen, 
experiment  1 15 


Uncorrected  and  corrected  stress  waves  for  a  12.7mm  specimen, 
experiment  134 


Figure  4. 


CHAPTER  5 

EXPERIMENTAL  RESULTS 


Introduction 

In  this  research  fifty  experiments  were  conducted  on  compacted  sand 
specimens.  The  compaction  moisture  and  density  combinations  ranged  from 
dry  of  optimum  to  wet  of  optimum  conditions  as  determined  by  the  Harvard 
miniature  compaction  procedure  (see  figure  3.2)  The  majority  of 
experiments  were  performed  with  specimens  prepared  at  or  near  optimum 
conditions.  A  minimum  of  two  experiments  were  conducted  at  each 

combination  of  specimen  moisture  content  and  density.  The 

moisture/ density  data  points  for  each  group  and  their  relationship  to 
the  Harvard  miniature  compaction  curve  are  presented  in  appendix  H.  The 
applied  stresses  and  maximum  strain-rates  for  the  experiments  ranged 

from  130  MPa  to  760  MPa,  and  from  500  s-1  to  5000  s-1  5  respectively. 

The  applied  stress,  peak  stress,  strain  at  peak  stress,  and  strain-rate 
at  peak  stress  for  each  experiment  are  tabulated  in  appendix  E. 

Three  different  length  striker  bars  were  available  to  generate  the 
applied  stress  wave,  0.127,  0.25 and  0.508  m.  Dispersion  effects 
could  not  be  resolved  for  the  0.508  m  bar  and  the  0.127  m  bar  would  not 
allow  the  specimen  to  reach  equilibrium  within  the  duration  of  the 
stress  wave  (»  50  microseconds).  Hence,  the  0.25*1  m  striker  was 
selected  to  initiate  the  stress  wa,Te  for  all  the  experiments  reported 
herein. 

Although  50  experiments  were  performed,  only  the  results  of  27  are 
used  for  the  purpose  of  comparison.  In  the  process  of  preparing  the 
SHPB  apparatus  for  firing,  a  seating  strain  was  applied  to  the 
specimen.  The  amount  of  seating  strain  sustained  by  a  specimen  varied 


from  experiment  to  experiment.  This  was  an  artifact  of  the  SHPB 
apparatus  and  required  a  change  in  the  experimental  procedure  so  that 
the  inherent  variability  was  accounted  for  in  the  data  reduction 
procedure.  The  cause  of  this  difficulty  and  the  procedure  that  was 
developed  to  measure  the  variable  seating  is  discussed  in  chapter  6. 

The  procedure  to  measure  the  seating  strain  was  applied  to 
experiment  112  and  every  experiment  thereafter.  Accurate  seating 
strains  are  necessary  to  determine  the  change  in  specimen  length  and 
phase  relationships  (e.g.,  initial  gas  porosity).  A  knowledge  of  these 
parameters  is  essential  as  they  will  influence  specimen  response  to  a 
given  stress.  If  accurate  comparisons  of  tests  are  to  be  made  the 
initial  length  of  the  specimens  being  compared  must  be  known.  Hence,  in 
the  experiments  where  seating  strains  were  not  measured  (i.e.,  31 
through  76),  the  results  are  not  suitable  for  the  purpose  of  comparison. 

The  average  seating  strain  experienced  by  the  specimens  in 
experiment  112  through  1 67  was  7.5  percent  with  a  standard  deviation  of 
3.8.  The  results  are  tabulated  in  table  E.4  of  appendix  E.  For  the 
experiments  where  the  seating  strain  was  not  measured,  the  initial 


specimen  length 

was 

adjusted  by 

7.5 

percent  for 

use  in 

computations. 

Figure 

5.1  shows 

the 

stress-strain 

response 

experiment  55  using  the  specimen  length  that  was  measured  just  before  it 
was  placed  between  the  bars  and  the  adjusted  length.  It  is  evident  that 
the  overall  response  trend  is  not  affected  by  adjusting  the  specimen 
length  by  7.5  percent.  While  the  true  response  to  use  for  comparison  is 
not  known,  the  response  computed  for  7.5  percent  is  not  inconsistent 
with  the  response  obtained  in  experiments  where  the  seating  strain  was 
measured.  Hence,  although  the  specimen  stress-strain  response  for 


experiments  31  through  76  is  consistent,  they  are  not  suitable  to 
validate  experimental  replication.  Plots  of  the  computed  stress-strain 
response  and  summary  information  for  each  of  the  50  experiments  are 


presented  in  appendix  I. 

Stress-Strain  Response 

Figure  5.2  shows  the  stress-strain  response  for  a  12.7  mm  specimen 
to  an  applied  stress  of  400  MPa.  This  typical  stress-strain  curve  can 
be  divided  into  three  regions  which  describe  the  overall  specimen 
response.  The  first,  region  0  to  A,  inertia  effects  act  to  oppose  the 
equilibration  of  stress.  The  extent  of  this  region  can  be  estimated  by 
using  the  criterion  of  Davies  and  Hunter  (1963)  as  well  as 
experimentally.  At  A  the  specimen  is  considered  to  be  in  quasi¬ 
equilibrium  and  the  deformation  nearly  homogenous.  Region  A  to  B  is 
characterized  by  the  closing  of  air  voids  and  rearrangement  of  particles 
into  a  denser  packing.  The  strain  at  B  is  equal  to  the  initial  gas 
porosity  of  the  specimen.  Computation  of  volumetric  phase  relationships 
indicates  that  at  this  point  the  air  volume  is  essentially  zero.  It  has 
been  shown  that  in  dynamic  uniaxial  strain  experiments  on  partially 
saturated  sand  where  the  stress  changes  are  large,  the  pore  air  will  be 
highly  compressed,  and  then  dissolved  in  the  pore  water  (Whitman, 
1970).  This  behavior  is  characteristic  of  the  response  observed  in 
region  A  to  B.  With  no  pore  air  remaining  (zero  gas  porosity),  the 
specimen  will  be  fully  saturated  (point  B).  Under  these  conditions  the 
compressive  resistance  of  the  water  will  greatly  exceed  the  compressive 
resistance  of  the  soil  skeleton  (Whitman,  1970)  and  the  specimen  will 
become  substantially  stiffer.  This  behavior  governs  the  response  in  the 
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region  beyond  B  up  to  nearly  the  peak  at  C. 

In  the  unloading  region  of  the  curve  the  specimen  in  many  of  the 
experiments  continued  to  acculumate  strain  (e.g.,  see  figures  4.13  and 
4.14).  The  cause  of  the  additional  strain  accumulation  may  be  that  the 
specimen  can  retain  some  memory  of  its  loading.  Another  cause  might  be 
the  extrusion  of  material.  The  potential  for  loss  of  material  during 
the  experiment  is  discussed  in  chapter  6.  Due  to  these  uncertainties, 
the  unloading  region  of  the  curve  has  not  been  considered  in  the 
analyses  reported  herein. 

Uniaxial  Strain  Path 


A  nearly  uniaxial  strain  state  is  forced  upon  the  specimen  by  the 
fact  that  it  is  contained  in  a  confining  cylinder.  This  condition  may 
be  verified  by  comparing  the  radial  strain  experienced  by  the  confining 
cylinder  and  maximum  longitudinal  strain  experienced  by  the  specimen. 
The  radial  strain  was  estimated  by  computing  a  maximum  radial  deflection 
of  the  confining  cylinder  using  thick-walled  cylinder  theory.  The 
radial  deflection  of  a  thick-walled  cylinder  (AR)  subject  to  internal 
pressure  may  be  found  by  (Seely  and  Smith,  1952); 


AR 


PR,  R?  ♦  R? 

~E  tR2  T  R2  +  V]  * 


(5.1) 


where  R1  is  the  inside  radius  of  the  cylinder,  R2  is  the  outside  radius 
of  the  cylinder,  P  is  the  internal  pressure,  E  (10  GPa)  is  Young's 
modulus  of  the  confining  cylinder,  and  v  (0.3)  is  Poisson's  ratio  for 


the  confining  cylinder.  The  internal  pressure  was  taken  as  the  peak 
stress  experienced  by  the  specimen  which  is  equivalent  to  an  assumption 


that  Poisson's  ratio  for  the  soil  is  0.5.  This  assumption  is  reasonable 


if  the  soil  is  near  saturation  as  the  pore  water  will  be  highly 
resistant  to  volume  change  (Whitman,  1970).  In  nearly  all  the 
experiments  conducted  volumetric  computations  indicated  that  the 
specimens  were  completely  saturated  at  peak  stress.  The  computed  radial 
strain  is  thought  to  be  an  upper  bound  as  the  internal  pressure  is  not 
distributed  over  the  entire  length  of  the  confining  cylinder  due  to  some 
overlap  on  the  pressure  bars  (see  figure  3.3).  The  computed  radial 
strains  are  presented  in  table  E.5  in  appendix  E. 

The  computed  radial  strain  of  the  confining  cylinder  for  all 
experiments  did  not  exceed  1  percent.  Comparing  the  radial  strain  of 
the  confining  cylinder  (er)  to  the  maximum  longitudinal 
strain  (e^)  experienced  by  the  specimen  (see  table  E.3  in  appendix  E) , 
it  was  found  that; 


e  <  0.06  c.  ,  (5.2) 
r  l 

for  all  experiments  except  one  (experiment  137),  where  the  radial  strain 
was  9  percent  of  the  longitudinal  strain.  Hence,  by  comparison,  the 
specimen  can  be  considered  to  be  in  a  state  of  nearly  uniaxial  strain 
during  the  experiment. 


Experimental  Replication 

If  the  SHPB  technique  is  to  prove  useful  for  evaluating  dynamic 
soil  response,  the  reproducibility  of  experimental  results  must  be 
established.  This  section  will  present  the  results  of  replicate 
experiments  on  different  length  specimens  subjected  to  a  range  of 


applied  stresses. 

Replicate  experiments  were  conducted  at  nominal  applied  stresses  of 
250,  400,  and  520  MPa  and  specimen  lengths  of  6.35,  and  12.7  mm.  Due  to 
variability  in  test  conditions  actual  stresses  varied  by  less  than  ±  8 
percent  from  the  nominal  stress.  The  results  are  presented  in  figures 
5.3  through  5.8  with  the  actual  applied  stress  given  in  parenthesis  next 

to  the  experiment  indentifi cation.  These  figures  show  that  the 

experimental  results  can  be  considered  to  be  reproducible.  It  can  be 
seen  that  in  nearly  all  cases  the  slopes  of  the  stress-strain  curves, 
the  peak  stress,  and  the  strain  at  peak  stress  for  the  replicated 
experiments  are  directly  comparable. 

However,  there  are  several  discrepancies  which  may  be  noted. 
Figure  5.5  for  6.35  mm  specimens  with  250  MPa  applied  stress  shows  both 

specimens  behaving  similarily  up  to  275  MPa.  At  that  point,  the 

specimen  in  experiment  1 62  began  to  accumulate  strain  with  very  small 
changes  in  stress,  whereas  the  specimen  in  experiment  163  continued  to 
accumulate  stress  with  very  small  changes  in  strain.  The  initial 
specimen  lengths,  phase  relationships,  and  applied  stresses  were 
essentially  identical  for  both  specimens.  Figure  5.7  for  12.7  mm 
specimens  with  400  MPa  applied  stress  shows  the  stress-strain  response 
for  each  specimen  to  be  very  similar  except  that  each  begins  to  stiffen 
at  different  strains. 

These  discrepancies  may  be  attributed  to  several  factors.  It  is 
difficult  to  prepare  specimens  with  indentical  properties  (i.e., 
moisture,  density,  etc).  Since  seating  strains  cannot  be  controlled 
with  great  precision  it  is  difficult  to  have  specimens  with  identical 
phase  relationships.  Variable  friction  force  between  the  launch  tube 


and  the  striker  bar  also  makes  it  difficult  to  achieve  duplicate  impact 
velocities.  These  potential  sources  of  difficulty  in  replicating 
experimental  results  are  discussed  in  chapter  6.  Another  factor  which 
might  hamper  the  replication  of  experimental  results  is  that  soil  and  or 
moisture  might  be  escaping  through  the  annular  region  between  the 
confining  cylinder  and  the  pressure  bars  during  the  experiment  (see 
figure  3.3)  A  discussion  of  soil  and  moisture  loss  is  presented  in 
chapter  6. 

Soil  Behavior 

The  observed  soil  behavior  may  be  characterized:  (1)  by  comparison 
of  the  stress-strain  response  of  specimens  with  similar  lengths  to  a 
range  of  applied  stresses,  (2)  by  comparing  the  response  of  specimens 
with  different  lengths  to  the  same  applied  stress,  and  (3)  by  comparing 
the  response  of  the  specimens  with  similar  lengths  but  prepared  at 
different  initial  moisture  contents  to  the  same  applied  stress. 

Figures  5.9  and  5.10  show  the  typical  stress-strain  response  for 
the  nominal  specimen  lengths  of  6.35  mm  and  12.7  mm  respectively,  to  a 
range  of  applied  stresses.  The  average  stress  experienced  by  each 
specimen  increased  with  increasing  applied  stress  independent  of  the 
specimen  length.  For  both  specimen  lengths,  the  stress-strain  response 
is  very  similar  for  applied  stresses  up  to  400  MPa  with  some  increase  in 
stiffness  observed  at  the  higher  (520  MPa)  applied  stress.  For  all 
applied  stress  levels,  the  specimens  began  to  stiffen  at  strains 
approximately  equal  to  the  initial  gas  porosity.  For  both  specimen 
lengths  and  at  all  applied  stresses,  the  strain  at  peak  stress 
experienced  by  the  specimen  exceeded  the  Initial  gas  porosity  of  the 


material. 


Figures  5.11,  5.12,  and  5.13  present  the  stress-strain  response  of 
specimens  to  the  same  approximate  applied  stress  but  with  different 
specimen  lengths.  The  nominal  applied  stresses  are  250,  400,  and  520 
MPa  respectively.  It  can  be  seen  (in  figures  5.11  and  5.12)  that  the 
strain  experienced  by  the  shorter  specimens  exceeded  those  of  the  longer 
specimens.  Also  the  longer  specimens  showed  greater  stiffness  at  lower 
strains  than  the  shorter  specimens. 

The  effect  of  moisture  content  variations  on  stress-strain 
response,  was  investigated  with  specimens  that  were  prepared  at  nominal 
moisture  contents  of  7,  13,  and  15  percent.  Figure  5.14  shows  the 
effect  of  moisture  content  on  the  stress-strain  response  of  6.35  mm 
specimens  to  a  nominal  applied  stress  of  250  MPa,  and  figure  5.15  shows 
the  response  of  12.7  mm  specimens  to  a  nominal  applied  stress  of  400 
MPa.  These  figures  show  that  the  average  stress  experienced  by  the 
specimens  increased  while  the  strain  at  peak  stress  decreased  with 
increasing  moisture  content.  As  with  the  other  stress-strain  curves 
shown,  there  is  a  marked  break  in  slope  near  a  strain  equal  to  the 
initial  gas  porosity.  This  change  in  slope  is  not  observed  for  the  12.7 
mm  specimen  with  the  lowest  moisture  content  (experiment  135)  because 
the  maximum  strain  (16  percent)  did  not  approach  the  initial  gas 
porosity  of  23.4  percent.  Also,  the  specimens  became  somewhat  stiffer 
with  increasing  moisture  content,  at  least  at  strains  in  excess  of  the 
initial  gas  porosity.  At  strains  less  than  the  initial  gas  porosity 
this  observation  is  less  clear. 
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Figure  5  &  Replicate  experiments  for  12.7mm  specimens  at  an  applied 
stress  of  approximately  250  MPa 


Figure  59  Stress  vs  strain  response  for  6  35mm  specimens 
of  applied  stresses 


145  (269  MPa) 


Figure  5  10  Stress  "s.  strain  response  for  12.7mm  specimens  to  a  range 
of  applied  stresses 


Figure  5  1 1  Comparison  of  stress  vs.  strain  response  based  on  specimen 
length  to  an  applied  stress  of  approximately  2  50  MPa 


Comparison  of  stress  vs.  strain  response  for  6.35mm  specimens 
based  on  moisture  content  to  an  applied  stress  of 
approximately  250  MPa 
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CHAPTER  6 


DISCUSSION  OF  EXPERIMENTAL  RESULTS 


Introduction 

Chapter  4  presented  the  method  used  to  develop  meaningful  stress- 
strain  data  from  the  raw  SHPB  data.  Chapter  5  presented  the  stress- 
strain  response  for  sand  specimens  compacted  to  different 
moisture/ density  conditions  and  subjected  to  a  range  of  applied 
stresses.  This  chapter  discusses  the  results  from  the  SHPB  tests  on  the 
sand  specimens  and  their  limitations.  It  will  also  discuss  the 
application  off  the  SHPB  technique  to  measure  the  high  strain-rate 
response  of  soils. 

Experimental  Assumptions 

From  the  strain  gauge  data  collected  during  a  SHPB  experiment  it  is 
assumed  that  the  stress,  strain,  and  strain-rate  response  in  the 
specimen  at  any  time  can  be  computed  by  averaging  the  forces  and 
particle  velocities  at  the  specimen-bar  interfaces.  For  the 
computations  to  be  meaningful  the  assumptions  of  the  experiment  must  be 
satisfied  (see  appendix  D).  The  degree  to  which  the  assumptions  were 
satisfied  in  the  tests  reported  herein  for  soil  specimens  and  large 
diameter  pressure  bars  is  also  presented  in  appendix  D. 

Experimental  Replication 

For  the  SHPB  technique  to  be  effective  in  measuring  soil 
properties,  its  results  must  be  reproducible.  Replicate  experiments 
were  conducted  on  essentially  equivalent  test  specimens  and  test 


conditions.  In  this  research,  several  aspects  of  the  experimental 
process  have  been  identified  as  hindering  complete  experimental 

reproducibility  with  compacted  sand  specimens. 

In  attempting  to  achieve  replicate  experimental  results,  each 
specimen  was  prepared  according  to  the  same  procedure.  The  procedure 
for  preparing  each  specimen  has  been  presented  in  chapter  3.  However, 
regardless  of  following  a  consistent  procedure  for  preparation, 
variations  in  specimen  parameters  existed  from  specimen  to  specimen. 

For  example,  the  optimum  moisture  content  for  the  soil  was  13.3  percent 
the  moisture  content  for  the  specimens  prepared  to  be  at  optimum  ranged 
from  10.4  to  13*9  percent.  The  results  shown  in  chapter  5  have 

demonstrated  that  the  initial  parameters  will  govern  the  specimen 

stress-strain  response  and  slight  variations  in  initial  parameters  have 
led  to  significant  difference  in  measured  properties. 

Due  to  the  nature  of  the  specific  SHPB  used  for  this  research  and 
the  compressibility  of  the  specimens,  the  seating  strains  produced  by 
positioning  the  specimen  between  the  bars  could  not  be  controlled. 
After  the  specimen  was  positioned  between  the  bars,  and  a  tight  contact 
between  system  components  checked,  the  hydraulic  actuator  was  used  to 
position  the  momentum  trap.  During  this  process,  a  variable  stress  was 
applied  through  the  transmitter  bar  to  the  specimen  causing  it  to 
compress.  The  range  of  this  stress  was  from  2  to  8  MPa.  A  measure  of 
the  seating  strain  was  obtained  by  recording  the  displacement  the 
momentum  trap  caused  by  the  movement  of  the  hydraulic  actuator.  However 
the  inability  to  control  the  amount  of  seating  strain  made  it  difficult 
to  have  two  specimens  with  identical  lengths  and  phase  relationships 
(e.g.,  initial  gas  porosity).  This  increased  the  difficulty  in 


achieving  complete  experimental  replication  and  contributed  to  the 
observed  discrepancies. 

Replication  of  experimental  results  is  also  dependent  on  how 
consistently  the  load  is  applied.  As  noted  in  chapter  3  the  magnitude 
of  the  applied  load  is  proportional  to  the  impact  velocity  of  the 
striker  bar  which  is  controlled  by  the  breech  pressure  in  the  gas  gun. 
Over  the  period  of  time  the  experiments  were  conducted,  the  breech 
pressure  to  produce  a  given  impact  velocity  differed  from  that  which  was 
predicted  from  the  theoretical  curve  due  to  a  variable  friction  force 
between  the  launch  tube  and  the  striker  bar.  To  minimize  the  effects  of 
friction,  the  launch  tube  was  periodically  swabbed  out  and  treated  with 
a  light  coating  of  lubricating  oil  and  a  dust  cover  was  fabricated  for 
the  vent  holes  on  the  launch  tube  so  that  contamination  could  be  reduced 
when  the  system  was  inactive.  Regardless  of  these  efforts,  it  was 
difficult  to  achieve  duplicate  impact  velocities  with  a  given  breech 
pressure. 

In  spite  of  the  above  mentioned  difficulties,  experimental 
replication  can  be  achieved  when  using  soil  specimens  in  a  SHPB 
experiment.  This  has  been  demonstrated  convincingly  by  figures  5.3,  5.6 
and  5.8. 

Soil  Behavior 

In  one-dimensional  compression  the  general  stress-strain  response 
exhibited  by  soil  is  S-shaped.  For  small  stresses  changes,  yielding  is 
observed  with  the  stress-strain  curve  concave  to  the  strain  axis.  For 
large  stress  changes,  the  behavior  is  character! zed  by  stiffening  with 
the  stress-strain  curve  being  concave  to  the  stress  axis.  This  behavior 


is  shown  in  figure  6.1.  The  soil  used  for  the  test  reported  in  this 
figure  is  a  dry  desert  alluvium  taken  from  a  site  near  Yuma,  Arizona 
(Dass  and  Bratton,  1983).  The  general  specimen  stress-strain  response 
observed  in  this  research  is  consistent  with  this  description  and  is 
similar  to  that  found  by  other  investigators  who  performed  uniaxial 
strain  experiments  on  similar  soils  but  a  lower  loading  rates  (Jackson, 
1968;  Calhoun  and  Kraft,  1966). 

In  figures  5.9  and  5.10,  at  strains  less  than  the  initial  gas 
porosity  the  compressibility  of  the  specimens  is  nearly  constant.  In 
figure  5.9  the  stress-strain  response  for  the  specimen  of  experiment  1 63 
does  not  adhere  to  this  behavior.  However,  the  initial  gas  porosity  of 
the  specimen  was  50  percent  less  than  that  of  the  other  two  so  that  a 
stiffer  response  at  a  lower  strain  is  to  be  expected.  At  strains  in 
excess  of  the  initial  gas  porosity,  the  compressibility  of  the  specimens 
is  again  nearly  constant  but  with  a  value  greater  than  the  initial 
compressibility.  For  example  in  figure  5.10  the  tangent  modulus  for 
experiment  13^  at  4  percent  strain  is  0.2  GPa  whereas  at  11  percent  it 
is  1.5  GPa.  It  should  be  noted  that  the  tangent  modulus  at  11  percent 
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compressibility  lies  a  transition  zone  which  can  be  identified  through 
the  initial  gas  porosity.  In  this  zone  the  response  of  the  specimen  is 
shifted  from  being  governed  by  the  soil  mass  to  being  governed  by  the 
pore  water.  Hence,  the  initial  gas  porosity  of  a  specimen  can  be  used 
as  an  indicator  to  determine  a  change  in  the  response  governing 
mechanism.  It  is  also  of  interest  to  note  that  a  similar  response  has 
been  observed  for  McCormick  Ranch  sand  subjected  to  high  hydrostatic 
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compression  (Mazanti  and  Holland,  1970). 

Figures  5.11,  5.12,  and  5.13  compared  the  response  of  different 
length  specimens  to  a  range  of  applied  stresses.  In  chapter  5  the 
discrepancies  between  the  response  for  the  two  specimen  lengths  were 
noted  and  an  examination  of  the  specimens  initial  gas  porosities  were 
proposed  as  a  means  to  explain  them.  In  figures  5.11  and  5.12  the 
initial  gas  porosity  of  the  longer  specimens  were  approximately  50 
percent  less  than  that  of  the  shorter  specimens  (i.e.,  6.2  percent  as 
compared  to  9.0  percent;  respectively  for  figure  5.10  and  5.8  percent  as 
compared  to  10.6  percent  respectively;  for  figure  5.11).  Hence,  greater 
stiffness  at  lower  strains  is  to  be  expected  for  the  longer  specimens. 
In  addition,  smaller  strains  should  be  expected  for  specimens  with  lower 
initial  gas  porosities.  In  figure  5.13  the  situation  is  reversed  with 
the  initial  gas  porosity  of  the  shorter  specimen  less  than  of  the  longer 

(i.e.,  ^.8  percent  as  compared  to  7.5  percent;  respectively).  Hence, 

the  behavior  shown  in  figure  5.13  is  consistent  with  the  conclusion  that 
the  specimen  stress-strain  response  is  governed  by  the  initial  gas 
porosity.  In  addition,  it  can  be  concluded  that  the  stress-strain 
response  of  specimens  with  different  lengths  to  the  same  applied  stress 
will  be  similar  if  their  initial  gas  porosities  are  similar. 

Figures  5. 11)  and  5.15  compare  the  stress-strair.  response  of 
specimens  with  similar  lengths  but  prepared  at  different  moisture 
contents  to  the  same  nominal  applied  stress.  The  response  shown  in 

these  figures  for  both  specimen  lengths  is  consistent  with  the  exception 

of  the  specimen  of  experiment  167.  However,  as  was  the  case  for 
experiment  163  discussed  above,  the  initial  gas  porosity  of  the  specimen 
is  approximately  50  percent  lower  than  the  other  two  specimens.  Hence 
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the  response  is  not  inconsistent  with  what  has  already  been  presented. 
For  the  specimens  prepared  near  and  wet  of  optimum,  the  stress-strain 
response  is  consistent  with  the  results  presented  throughout  chapter  5; 
hence  their  analysis  will  follow  what  has  already  been  discussed. 

In  figures  5. 11*  and  5.15  the  experiments  where  the  specimens  were 
prepared  dry  of  optimum  were  135  and  139.  For  the  specimen  of 
experiment  135  a  change  in  slope  was  not  observed  and  the  strain  did  not 
approach  the  initial  gas  porosity.  It  can  also  be  seen  that  the  slope 
of  the  stress-strain  curve  is  nearly  constant  which  indicates  that  the 
compressibility  of  the  soil  skeleton  can  be  considered  to  be  nearly 
constant  and  that  the  compression  of  the  air  voids  is  directly 
proportional  to  changes  in  stress. 

A  difference  is  observed  in  the  response  of  the  specimens  prepared 
at  the  moisture/density  conditions  dry  of  optimum  for  the  two  specimen 
lengths  In  the  12.7  mm  specimen  (experiment  135),  which  had  a  gas 
porosity  of  23- 4  percent,  very  little  build  up  of  stress  was  experienced 
at  the  maximum  strain  of  16  percent.  For  the  6.35  mm  specimen 
(experiment  139),  which  had  a  gas  porosity  of  17.7  percent,  a 
substantially  higher  stress  was  built  up  as  the  strain  approached  this 
value.  Two  factors  may  be  responsible  for  this  difference  in 
behavior.  First,  the  density  of  the  6.35  mm  specimen  (1.98  g/cm^)  is 
greater  than  that  of  the  12.7  mm  specimen  (1.87  g/cm^).  Second,  the 
initial  gas  porosity  of  the  6.35  mm  specimen  is  lower  than  that  of  the 
12.7  mm  specimen.  As  the  effect  of  density  variations  on  specimen 
response  was  not  investigated  it  can  not  be  concluded  that  one  factor  is 
more  dominant  than  the  other.  However,  the  response  of  the  6.35  mm 
specimen  is  consistent  with  the  response  of  other  specimens  where  the 
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strain  approached  or  slightly  exceeded  the  initial  gas  porosity  (e.g., 
experiment  119);  hence  supporting  the  conclusion  that  the  initial  gas 
porosity  is  a  governing  parameter  in  determining  specimen  stress-strain 
response. 

Strain  in  Excess  of  Initial  Gas  Porosity 

As  discussed  in  section  5.2  the  initial  gas  porosity  is  an 
approximate  upper  bound  on  the  strain  that  a  specimen  can  sustain. 
However,  the  results  presented  in  section  5.5  showed  that  for  nearly 
every  experiment,  the  strain  sustained  by  the  specimen  exceeded  the 
initial  gas  porosity  (see  figure  6.2).  There  are  several  factors  which 
may  be  working  together  to  account  for  this  discrepency: 

(1)  loss  of  soil  and  or  moisture, 

(2)  compression  of  the  pore  water  as  the  specimen  becomes 
saturated,  and 

(3)  radial  expansion  of  the  confining  cylinder  and  specimen. 

These  factors  will  be  examined  in  the  following  paragraphs. 

Soil  and  moisture  loss  are  possible  because  the  diameter  of  the 
bars  and  the  inside  diameter  of  the  confining  cylinder  are  not  equal. 
Although  the  tolerance  is  very  small  (see  figure  3*3)*  soil  and  moisture 
loss  can  occur  through  the  annular  region  between  them.  The  soil  loss 
is  induced  by  a  pressure  gradient  which  is  set  up  between  specimen-bar 
interfaces  and  the  atmosphere. 

At  the  specimen- incident  bar  interface  a  portion  of  the  initial 
compressive  stress  wave  is  reflected  as  a  tensile  wave  due  to  the  lower 
impedance  of  the  specimen  relative  to  the  bar.  This  tensile  wave 
travels  back  down  the  incident  bar  toward  the  end  at  which  the  Impact 


occur ed.  As  the  impact  end  of  the  incident  bar  is  now  a  free  end,  the 
tensile  wave  is  reflected  as  a  compressive  wave  travelling  once  again 
toward  the  specimen,  hence,  reloading  the  specimen.  For  the  SHPB 
experiment  the  specimen  response  of  interest  is  that  due  only  to  the 
initial  stress  wave. 

The  presence  of  multiple  impacts  complicates  the  ability  to 

quantify  the  amount  of  soil  and  moisture  lost,  as  the  quantity  lost 
during  each  Impact  is  indeterminate.  In  addition,  the  specimen  may 
sustain  additional  strain  during  the  second  impact,  making  it 
meaningless  to  compare  a  measured  strain  with  the  computed  strain  which 
results  from  the  initial  stress  wave. 

An  attempt  to  measure  the  mass  of  soil  lost  was  made  in  several 
experiments.  The  mass  of  soil  lost  was  determined  by  recording  the 

combined  mass  of  the  confining  cylinder  and  specimen  before  it  was 

positioned  between  the  bars  and  just  after  it  was  removed,  at  the 
completion  of  the  experiment,  and  then  taking  the  difference.  The 

maximum  and  minimum  soil  loss  was  computed  to  be  7.9  g  and  1.7  g 

respectively.  The  results  are  tabulated  in  table  E.6  of  appendix  E. 
These  measurements  overestimate  the  amount  of  soil  lost  during  the 

analyzed  impact  because  the  specimen  is  subjected  to  multiple  impacts. 

Because  a  moisture  content  was  determined  at  the  time  the  specimen 
was  prepared  and  immediately  after  the  experiment,  a  measure  of  the 

amount  of  moisture  loss  sustained  by  the  specimen  can  be  computed.  The 

maximum  and  minimum  moisture  losses  were  computed  to  be  27.6  and  1.7 
percent  of  the  initial  moisture  content;  respectively.  The  moisture 
loss  for  each  experiment  is  tabulated  in  table  E.7  of  appendix  E.  If 
all  the  moisture  loss  is  assumed  to  occur  during  the  first  impact,  the 


specimen  can  sustain  an  additional  strain  proportional  to  the  volume  of 
moisture  lost.  Using  this  assumption,  the  strain  contribution  from 
moisture  loss  for  each  experiment  has  been  computed  and  the  results 
tabulated  in  table  6.1.  These  strains  will  also  be  overestimates  since 
multiple  impacts  occur  in  every  test. 

In  most  engineering  applications,  water  is  considered  to  be  an 
incompressible  fluid.  At  the  stresses  applied  in  this  research  (>  100 
MPa),  this  assumption  is  not  valid.  The  compression  of  the  pore  water 
will  contribute  to  the  strain  sustained  by  the  specimen. 

The  specific  volume  of  water  at  a  variety  of  temperatures  subjected 
to  a  range  of  applied  stresses  has  been  tabulated  in  the  Handbook  of 
Physical  Constants  (Sydney,  1966).  These  values  have  been  used  to 
compute  the  strain  contribution  due  to  the  compression  of  the  pore 
water.  The  computation  results  are  presented  in  table  6.1.  The  stress 
used  in  the  computations  was  taken  as  the  peak  stress  sustained  by  the 
specimen. 

It  has  been  shown  in  chapter  5  that  the  specimens  are  constrained 
to  a  state  of  nearly  uniaxial  strain  because  of  the  small  amount  of 
radial  strain  that  occurs.  However  small,  the  radial  strain  will  make 
some  contribution  to  the  overall  specimen  strain.  The  strain 
contribution  is  computed  by  accounting  for  the  change  in  specimen  area 
resulting  from  the  radial  displacement. 

Table  6.1  sums  the  strain  contributions  from  water  compression, 
moisture  loss,  and  radial  strain  and  compares  it  to  the  specimen  strain 
at  peak  stress.  These  data  are  plotted  in  figure  6.3*  It  can  be  seen 
that  in  all  but  a  few  experiments  a  balance  of  strain  still  remains. 
Had  sufficient  information  been  obtained  to  compute  the  strain 


contribution  from  soil  loss  the  remaining  balance  for  each  experiment 
would  be  less. 


Due  to 

the  uncertainty 

in 

the  evaluation 

of  the 

strain 

contributions 

from  the  above 

mentioned  factors, 

a  gross 

strain 

adjustment,  to 

account  for  the  excess 

between  strain  at 

peak  stress,  and 

the  initial  gas  porosity  has 

been 

computed .  The 

details 

of  the 

computations  are  presented  in  appendix  J.  The  gross  strain  correction 
requires  none  of  the  assumptions  made  in  the  computations  of  the  strain 
contribution  from  moisture  loss,  pore  water  compression,  or  radial 
expansion.  It  is  based  solely  on  the  peak  stress  sustained  by  the 
specimen.  In  addition,  the  strain  contribution  of  soil  loss  is  al3o 
included. 

The  results  show  that  the  strain  balance  from  the  gross  correction 
is  less  than  that  computed  by  summing  the  individual  components 
indicating  that  the  strain  contribution  due  to  soil  loss  can  be 
significant  (see  figure  6.4).  Regardless  of  the  method  used,  an 
accounting  for  the  excess  between  the  strain  at  peak  stress  and  the 
initial  gas  porosity  can  be  made. 

Strain-Rate  Dependence 

To  determine  the  strain-rate  sensitivity  of  the  compacted  sand, 
stress-strain-rate  curves  at  constant  strains  were  constructed  for  each 
nominal  specimen  length.  The  data  were  taken  from  experiments  conducted 
at  the  nominal  applied  stresses  of  250  MPa  and  400  MPa.  The  results  are 
presented  in  figures  6.5  and  6.6  for  the  6.35  mm  and  the  12.7  mm  long 
specimens  compacted  at  moisture  and  density  conditions  near  optimum; 
respectively.  The  details  of  the  construction  are  as  follows. 
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For  each  experiment,  the  stress  and  strain-rate  for  a  given  strain 
was  recorded.  For  the  6.35  mm  specimens  the  range  of  strains  was  from  5 
percent  to  14  percent,  and  for  the  12.7  mm  specimens  the  range  of 
strains  was  from  5  percent  to  10  percent.  The  stress-strain-rate  data 
points  from  the  experiments  for  a  given  strain  are  then  grouped 
according  to  the  applied  stress  at  which  the  experiment  was  conducted 
and  then  plotted.  From  figures  6.5  and  6.6  it  can  be  seen  that 
considerable  scatter  results.  To  simplify  the  presentation,  the  group 
of  stress-strain-rate  data  points  for  a  given  strain  and  nominal  applied 
stress  were  averaged.  In  figures  6.5  and  6.6,  the  dashed  lines  are  the 
average  stress-strain-rate  trajectories  for  a  given  applied  stress.  The 
solid  lines  connect  points  of  constant  strain  between  the  stress-strain- 
rate  trajectories. 

In  this  type  of  plot,  the  degree  to  which  the  compacted  sand  is 

strain-rate  dependent  can  be  judged  by  the  slope  of  the  constant  strain 

curves  connecting  the  stress-strain-rate  trajectories.  If  the  slope  is 
zero,  it  can  be  concluded  that  the  compacted  sand  response  is  not 

strain-rate  dependent.  The  results  presented  show  that  the  response  of 
the  compacted  sand  is  not  dependent  on  strain-rate  for  strains  below  the 
initial  gas  porosity.  This  is  consistent  with  the  results  reported  by 
Gaffney  et  al.  (1985)  who  found  no  strain-rate  dependence  in  a  drier  but 
otherwise  similar  soil  for  total  strains  less  than  the  initial  gas 

porosity. 

It  was  observed  in  figures  6.11  and  6.12  that  the  stress  which 
produced  a  given  strain  did  not  increase  with  increased  strain-rate, 
indicating  that  the  specimen  response  was  not  dependent  on  strain-rate 
for  strains  less  than  the  initial  gas  porosity.  However,  this  apparent 


strain-rate  independence  should  be  viewed  with  caution  for  two  reasons 


first,  the  constant  strain  curves  used  to  show  this  apparent  strain-rate 
independence  are  only  rough  averages  developed  from  a  few  data  points, 
and  second,  the  factors  identified  to  account  for  the  discrepancy 
between  the  strain  at  peak  stress  and  initial  gas  porosity  cannot  be 
quantified  with  the  necessary  accuracy  to  determine  what  their  effect  on 
the  apparent  strain-rate  dependency  might  be. 

Conditions  of  Experiment 

In  order  to  satisfy  the  assumptions  of  the  experiment  (see  appendix 
D)  the  soil  specimens  were  prepared  with  aspect  ratios  less  than  or 
equal  to  0.2  and  contained  in  a  confining  cylinder.  The  small  aspect 
ratios  were  necessary  to  satisfy  the  uniform  stress  assumption  while  the 
containment  facilitated  the  control  of  boundary  conditions. 

If  the  specimens  had  been  longer,  stress  equilibrium  would  not  have 
been  achieved  during  the  early  portion  of  the  experiment.  Had  this  been 
the  case,  a  greater  portion  of  the  stress-strain  curve  would  have  had  to 
be  discarded  as  inertial  forces  would  have  domainated  the  response. 
Therefore,  soil  specimens  in  the  SHPB  experiments  are  restricted  to 
small  aspect  ratios.  It  may  be  possible  to  use  greater  aspect  ratios  if 
the  duration  of  the  applied  stress  wave  is  extended  (e.g.,  by  using  a 
longer  striker  bar). 

In  the  traditional  SHPB  experiment  the  specimen  is  in  a  state  of 
nearly  uniaxial  stress.  Because  the  unconflned  strength  of  soil  (e.g., 
<  0.1  MPa)  is  insignificant  compared  to  the  applied  stresses  (e.g.,  > 
100  MPa),  reproducing  the  traditional  SHPB  experiment  with  soil 
specimens  is  not  feasible.  Hence,  some  means  of  radial  confinement  is 


required.  Three  degrees  of  radial  support  can  be  considered:  minimal, 
intermediate,  and  stiff.  Providing  minimal  support  would  most  closely 
approximate  the  uniaxial  stress  state  in  the  traditional  experiment. 
However,  even  the  flimsiest  of  containment  systems  would  provide  a 
radial  restraint  that  would  be  significant  relative  to  the  unconfined 
Strength  of  the  soil.  If  an  intermediate  confinement  is  used,  radial 
stress  and  strain  histories  would  have  to  be  measured  to  evaluate 
specimen  response.  Therefore,  the  use  of  minimal  or  intermediate 
containment  systems  would  complicate  the  data  analysis.  In  contrast  to 
the  first  two,  a  stiff  confinement  system  is  easy  to  implement  and  has 
the  advantage  that  it  constrains  the  specimen  to  a  nearly  uniaxial 
strain  path  which  can  be  easily  duplicated  at  both  higher  and  lower 
strain-rates. 

Applications  of  the  SHPB  Technique 

In  recent  years  complex  constitutive  models  have  been  developed  for 
use  in  numerical  computations  to  predict  the  response  of  soil  in  the 
region  close  to  a  high  yield  weapon  explosion  (e.g.,  Rubin  and  Sandler, 
1977).  These  models  are  used  in  large  computer  codes  to  study  ground- 
shock  problems  of  both  the  free-field  and  soil-structure  interaction 
types  (Nelson,  Baron,  and  Sandler,  1971).  If  these  models  are  to 
accurately  predict  field  response  they  must  be  fit  to  data  which  reflect 
the  type  of  confinement,  magnitude  of  stress  change,  and  time  scale  of 
stress  application  to  be  expected  (Whitman,  1970).  These  requirements 
have  placed  demands  on  expermentalist  to  develop  laboratory  procedures 
that  will  reflect  field  conditions  and  provide  data  for  the  evaluation 
of  the  material  constants  which  appear  in  the  models  (Nelson  et  al., 
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1971) 


The  soil  test  most  often  selected  to  duplicate  the  loading 
condition  in  the  overpressure  region  associated  with  weapons  effects  is 
the  uniaxial  strain  test.  For  the  past  30  years  the  dynamic  uniaxial 
strain  device  has  been  used  to  investigate  the  fundamental  nature  of 
soil  behavior  to  rapidly  applied  loads.  Although  current  devices  can 
apply  stresses  up  to  400  MPa  with  loading  time  of  0.3  milliseconds  the 
need  for  a  laboratory  test  that  more  closely  simulates  the  environment 
near  a  weapon  explosion  is  clear.  The  SHPB  technique  is  a  means  to 
satisfy  the  demand  for  more  reliable  data  in  support  of  a  wide  variety 
of  weapon  effects  problems. 

It  is  common  practice  to  extrapolate  laboratory  data  to  evaluate 
the  material  constants  which  appear  in  a  constitutive  model  that  will  be 
used  to  represent  the  soil  behavior  to  an  explosive  event.  As  the  SHPB 
technique  can  be  used  to  evaluate  material  response  to  stresses  and 
strain-rates  that  are  closer  in  situ  conditions,  extrapolation  can  be 
avoided  and  the  material  constants  determined  directly.  This  can  lead 
to  a  decrease  in  the  time  required  to  evaluate  the  adequacy  of  a  model 
to  a  particular  problem  as  well  as  a  better  model  because  the  constants 
can  be  evaluated  from  data  that  more  closely  duplicates  field 
conditions.  In  addition,  a  reduction  in  field  costs  may  also  be 
realized  as  extra  in  situ  measurements  can  be  avoided  if  the  laboratory 
data  can  provide  a  closer  approximation  of  the  anticipated  field 
conditions . 

The  flyer  plate  experiment  is  a  particular  example  where  the  data 
from  a  SHPB  experiment  would  prove  useful.  The  purpose  the  flyer  plate 
experiment  is  to  develop  a  technique  for  field  calibrating  soil  stress 
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and  motion  instrumentation  in  the  pressure  range  of  100  to  1000  MPa.  In 
this  experiment  a  planar  impact  is  applied  to  an  instrumented  soil  test 
bed  by  an  explosively  driven  metal  plate  such  that  the  loading  path  is 
approximately  uniaxial  strain.  To  calibrate  the  instrumentation  the 
soil  must  be  properly  characterized  such  that  the  material  model  will 
accurately  predict  the  stress  and  motion  fields.  The  SHPB  technique  is 
the  most  appropriate  means  to  evaluate  the  required  material  properties 
as  precise  and  accurate  measurements  can  be  made  over  the  stress  range 
of  interest. 

The  SHPB  technique  can  also  be  used  to  extend  the  knowledge  of 
dynamic  soil  response  in  several  other  areas.  As  previously  described, 
during  the  experiment  the  specimen  was  subjected  to  multiple  impacts. 
This  phenomenon  may  be  exploited  to  investigate  soil  response  to 
multiple  loadings  in  rapid  succession  without  altering  the  specimen 
properties  in  between  loadings.  Hence,  providing  a  laboratory  means  to 
investigate  soil  response  to  a  sucessive  impulse  loading  environment. 
Another  area  of  interest  is  the  effect  of  confinment  on  soil  response. 
By  varying  the  degree  of  radial  support  to  a  soil  specimen  the  effect  of 
confinment  on  the  soil-strain  response  can  be  evaluated.  Such 
information  would  be  valuable  to  the  understanding  of  soil  response  ro 
airblast  loading  during  the  transition  from  the  superseismic  region  to 
the  outrunning  region.  Also  the  SHPB  technique  is  capable  of  providing 
information  on  stress  attenuation  and  changes  in  waveform  with  distance 
and  their  influence  on  stress-strain  response  as  a  function  of  stress 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusion 

The  analyses  contained  in  this  research  permit  the  following 
conclusions  regarding  the  use  of  the  SHPB  experimental  method  to  measure 
the  dynamic  response  of  soil: 

(1)  For  compacted  sand  specimens  constrained  to  a  nearly  uniaxial 
strain  state,  the  assumptions  necessary  to  obtain  meaningful 
data  from  a  SHPB  experiment  can  be  satisfied. 

(2)  Experimental  replication  can  be  achieved.  However,  the 
experimental  results  are  strongly  influenced  by  initial  weight- 
volume  relationships  in  the  test  specimens  (e.g.,  initial  gas 
porosity),  variations  induced  by  equipment  (e.g.,  seating 
stresses  and  barren  friction),  and  material  losses  (both  soil 
and  moisture)  during  testing. 

(3)  The  stress-strain  response  is  governed  principally  by  the 
initial  gas  porosity  of  the  specimen.  At  strains  less  than  the 
initial  gas  pororsity  the  specimen  compresses,  closing  the  air 
voids.  At  a  strain  near  the  initial  gas  porosity  the  specimen 
begins  to  stiffen  significantly. 

(4)  The  specimens  prepared  near  optimum  conditions  were  insensitive 
to  strain-rate  at  strains  less  than  the  initial  gas  porosity. 
At  strains  in  excess  of  the  initial  gas  porosity  a  strain-rate 
dependency  is  apparent.  Pore  water  pressure  is  likely  to  play 
a  role  in  this  change,  but  how  is  not  clear. 
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(5)  It  is  concluded  that  the  SHPB  experimental  method  can  be 
developed  to  determine  the  dynamic  response  of  soil.  Although 
at  present  a  few  limitations  are  associated  with  using  soil 
specimens  in  a  SHPB  experiment,  these  obstacles  appear  to  be 
capable  of  resolution.  The  SHPB  method  is  a  relatively  simple 
experiment  and  fully  developed  can  significantly  extend  the 
range  of  stresses  and  strain-rates  that  can  be  applied  beyond 
the  capabilities  of  current  equipment  used  for  dynamic  soil 
investigations. 

Recommendations 

The  success  of  this  research  effort  provides  a  justification  for 
the  expanded  use  of  the  SHPB  technique  for  evaluating  dynamic  soil 
response.  Consequently,  it  is  recommended  that  study  be  initiated  in 
three  areas:  (1)  alternative  measurements,  (2)  containment  systems,  and 
(3)  modification  to  equipment. 

As  a  supplement  to  the  traditional  measurement  techniques  which 
have  been  described  in  chapters  2  and  3,  the  feasibility  of  sensing  the 
stress  waves  directly  on  the  specimen-bar  interfaces  should  be 
explored.  For  foam  specimens  such  a  technique  has  been  employed  that 
uses  quartz  crystals  mounted  on  the  ends  of  the  pressure  bars  in  contact 
with  specimen  (Hodge  and  Wasley,  1969). 

The  soil  containment  system  used  in  this  research  has  presented 
several  difficulties  the  most  Important  of  which  is  that  loss  of  soil 
and  moisture  cannot  be  prevented.  Therefore  it  would  be  advantageous  to 
have  a  system  that  would  prevent  such  losses.  Also  such  a  system  would 
allow  unloading  behavior  to  be  investigated.  In  addition,  if  a 


confining  cell  similar  to  those  used  for  triaxial  testing  could  be  built 
that  would  be  able  to  apply  confining  stresses  of  several  hundred 
megapascals,  a  means  would  be  available  to  investigate  dynamic  shear 
response.  Such  a  containment  system  would  extend  the  capabilities  of 
the  method  when  using  soil  specimens  considerably.  As  a  minimum 
however,  measurements  of  radial  stresses  and  strain  should  be 
attempted.  This  information  alone  would  be  very  valuable  in  deriving 
constitutive  models  for  general  deformation  paths. 

The  final  area  in  which  recommendations  will  be  made  concerns 
experimental  equipment.  Two  main  problems  were  encountered  in  using  the 
equipment  described  in  chapter  3  when  using  soil  specimens:  (1)  the 
inability  to  accurately  control  seating  strains,  and  (2)  only  stresses 
in  excess  of  200  MPa  could  be  applied  with  any  degree  of  consistency. 

The  first  difficulty  can  be  eliminated  by  the  introduction  of 
controls  sensitive  to  small  changes  in  stress.  These  would  be 
beneficial  for  two  reasons.  First,  an  increased  control  on  seating 
strains  could  be  established  and  second,  a  preload  could  be  precisely 
controlled  to  simulate  geostatic  overburden  stresses. 

The  second  problem  may  be  overcome  by  using  a  striker  bar  with  a 
characteristic  impedance  less  than  that  of  the  steel  pressure  bars.  As 
the  applied  stress  is  proportional  to  the  characteristic  impedance  of 
the  striker  bar  (see  appendix  C)  a  reduction  in  the  impedance  will  cause 
a  smaller  3tress  to  be  applied  for  the  same  impact  velocity.  For 
example,  if  the  striker  bar  is  constructed  of  aluminum  which  has  a 
characteristic  impedance  nearly  three  times  less  that  that  of  steel,  for 
the  same  impact  velocity  the  applied  stress  would  be  reduced  by  nearly  a 


factor  of  three. 


A  final  recommendation  concerning  equipment  pertains  to  the 


pressure  bars.  The  ratio  of  the  characteristic  impedances  for  the  soil 
specimens  and  the  pressure  bars  used  in  this  research  was  approximately 
60  to  1 .  Although  the  uniform  stress  assumption  was  not  seriously 
compromised,  stress  equilibrium  in  the  specimen  would  have  been  achieved 
sooner  if  the  impedance  mismatch  had  been  smaller.  Also  the  constancy 
of  strain-rate  could  have  been  improved  if  the  characteristic  impedance 
of  the  pressure  bars  had  been  closer  to  that  of  the  soil  (Rand,  1967). 
Therefore  it  is  recommended  that  pressure  bars  be  constructed  of 
material  other  than  steel.  Potential  candidates  may  be  aluminum  and 
possibly  some  type  of  plastic. 

With  the  development  of  improved  equipment  and  testing  procedures, 
refinement  of  the  soil  behavior  measured  will  also  be  achieved. 
Specifically  the  nature  of  the  effects  of  repeated  impacts  and  their 
analysis  should  be  developed  as  well  as  the  characterization  of  the  soil 
behavior  of  the  soil  near  the  peak  of  the  stress-strain  curve.  The 
correlation  of  SHPB  data  with  more  conventional  dynamic  soil  tests  and 
field  experiments  should  also  be  pursued. 

While  the  work  reported  herein  has  established  the  potential  for 
application  ef  the  SHPB  to  critical  problems  in  dynamic  soil  behavior, 
further  research  and  development  as  outlined  above  will  be  necessary 
before  the  SHPB  results  may  be  used  routinely  or  with  any  substantial 
degree  of  confidence. 
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A.l  General 


The  procedure  for  correcting  a  waveform  for 
dispersion  can  be  broken  down  into  four  basic  operations 
(Follansbee  and  Frantz,  1983): 

(1)  transformation  of  the  original  waveform  to  a  Fourier 
series, 

(2)  computation  of  the  phase  velocity  for  each  frequency 
component, 

(3)  computation  of  the  adjustment  to  the  phase  angle 
based  on  position, 

(4)  reconstruction  of  the  original  waveform  with 
adjusted  phase  angles. 

A. 2  Transformation  to  a  Fourier  Series 

The  general  expression  for  a  Fourier  series  can  be 
written  as  follows  (Wiley  and  Barrett,  1982); 

oo 

f(t)  «  Ao  +  Z  An  cos  niUot  +  ®n  sin  n  coot  ,  (A.l) 

~2  n=i 

where  f (t)  is  a  periodic  function,  u>o  is  the  wave 
frequency  defined  as; 

U)0*  JJL  •  (A. 2) 

T 

where  T  is  the  period  of  the  wave,  and  Ac,  An,  and  Bn  are 
the  Fourier  coefficients  which  can  be  expressed  as; 


I 

> 


fT 

A  o  -  JL  /f  <  t>  dt  #  (A. 3) 

T  -»0 


An 

-  A  ] 

(rit, 

cos  n  oi0t 

dt. 

(A. 4) 

Bn 

T  J 

-  A  i 

0 

rT 

ff(t) 

sin  nu)0t 

dt. 

(A. 5) 

T  J 

0 

Upon  evaluation  of  the  Fourier  coefficients  A 0,  An,  and 
Bn*  the  general  Fourier  series  (equation  A.l)  is 
transformed  to  a  pure  Fourier  cosine  series.  The  reduction 
of  the  general  Fourier  series  involving  the  sum  of  a  sine 
and  cosine  term  is  accomplished  by  first  multiplying  the 
sum; 


A„  cos  nwot  +  B»  sin  n  u>0t  , 


(A. 6) 


by, 


V  An  +  Bn 

V  An  +  ®n 


to  yield. 


VAn  +  Bn 


An  cos  nuiot  +  Bn  sin  n  cu0t 


VAn  +  Bn 


VA  n  ♦  Bn 


By  establishing  the  new  variables; 


Dn  *  "VAn  ♦  ®n  • 


which  is  the  resultant  amplitude  of  the  frequency 


(A. 7) 


(A. 8) 


(A. 9) 
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component  nuiot«  end; 


cos<J>  ■  An  ,  (A.  10) 

Va„  +  B„ 

sin<t>  »  Bn  »  (A. 11) 

V An  +  Bn 

where  <J>  is  the  phase  angle  which  is  a  measure  of  the  lead 
or  lag  of  the  nth  harmonic  in  reference  to  the  sine  or 
cosine  wave  of  the  same  frequency  (Wiley  and  Barrett, 

1982).  Expression  (A. 8)  can  be  rewritten  as; 

Dn  (cos<|>cos  nuj0t  +  sin<J>sin  niu0t).  (A. 12) 

Applying  the  following  trigonometric  substitution  to 
expression  (A. 12)  (Selby,  1975); 

cos  a  cos  b  +  sin  a  sin  b  =  cos  (a  -  b)  ,  (A. 13) 

yields. 


Dn  cos  (nuuot  ”0)»  (A. 14) 

Substituting  expression  (A. 14)  into  equation  (A.l),  yields 
a  pure  Fourier  cosine  series  of  the  form; 

oo 

f(t)  «  Ao  +  Z  Dn  cos  (nu)0t  -0)  .  (A. 15) 

2  n=i 


There  is  one  note  on  the  computation  0 .  Instead  of 
computing  ip  as; 


<f>  vill  be  computed  as; 


<)>  *  cos  *  An  •  (A. 17) 

Dn 

This  is  because  the  range  of  the  arctangent  is  +.1T  /2, 
hence  for  some  phase  angles  (e.g.,  3  7T/4)  an  incorrect 
value  would  be  returned.  The  range  of  the  arccosine  is  0 
to  ir .  To  maintain  the  proper  sign  of  the  phase  angle  in 
the  computations;  if  Bn  <  0,  the  resultant  phase  angle 
will  be  the  inverse  sign  of  the  value  obtained  in  equation 
(A. 17) . 

A. 3  Computation  of  Phase  Velocity 

The  propagation  velocity  of  a  wave  at  a  given 
frequency  is  called  its  phase  velocity  (Cn).  Cn  can  be 
written  in  terms  of  cyclic  frequency  (f)  and  wavelength 
( A  )  as; 

Cn  ■  fnA  r  (A. 18) 

which  can  be  related  to  the  natural  frequency  (U>o) 
through; 


Equation  (A. 20)  can  be  written  in  dimensionless  form  by 
rearranging  and  dividing  both  sides  of  the  equation  by  C<>R 
to  yield. 


2T  Cn  R  »  R  n  U)o  .  (A. 21) 

A  Co  Co 

The  terms  on  the  right  side  of  equation  (A. 21)  are  all 
known  quantities.  The  only  unknown  in  equation  (A. 21)  is 
Cn*  A  value  for  Cn  can  be  computed  using  the  polynomial 
approximation  (equation  2.43  to  the  fundamental  mode  of 
the  dispersion  equation  (equation  2.42).  The  two  variables 
in  the  polynomial  approximation  are  Cn/Co ,  and  R/A  .  By 
computing  Cn/C0  for  a  given  value  of  R/A  ,  the  quantities 
on  each  side  of  equation  (A. 21)  can  be  compared.  Hence,  by 
iterating  through  values  of  R/A  ,  a  value  for  Cn/Co  can  be 
computed  such  that  the  two  sides  of  equation  (A. 21)  will 
be  equal.  Once  equality  has  been  established,  Cn  can  be 
computed  directly. 

A. 4  Phase  Angle  Adjustment  Based  On  Position 

From  the  dispersion  equation  (equation  2.42)  it  has 
been  determined  that  Cn  is  dependent  on  A  (Davies,  1948). 
This  causes  a  wave  to  disperse  as  it  propagates  because 
the  high  frequency  components  will  lag  behind  the  low 
frequency  components.  Follansbee  and  Frantz  (1983)  have 


derived  a  relationship  based  on  the  phase  angle  (0)  to 
account  for  the  dispersion  that  occurs  as  a  wave 
propagates  a  distance  Ax.  The  derivation  of  the  phase 
angle  correction  is  shown  below. 

The  term  inside  the  brackets  of  equation  (A. 15)  can 
be  rewritten  in  terms  of  position  as; 

nu)0t  -0  -  nu»0  Ft  +  Ax  -  A xl  ,  (A. 22) 

l  Co  CnJ 

where  Ax/Cn  is  the  time  for  a  given  component  to  travel  a 
distance  Ax.  Rewriting  the  terms  of  equation  (A. 22) 
yields ; 

n  Wq  t  -  0  -  nu*0 

nwgt  -  0  -  nu)g  t  -  ntun  A  X  f  Co  -  l"|  (A. 24) 

Co  l c„  J 

The  phase  angle  can  now  be  computed  as  a  function  of 
position  according  to; 

0/ ■  n  CUn  A  X  [  Cfl  -  ll  .  (A. 25) 

C0  [ cn  J 

A. 5  Wave  Reconstruction  „ 

The  last  operation  in  the  dispersion  correction 


t  + 


Ax 


(A. 23) 


procedure  is  wave  reconstruction.  Essentially  this  entails 
adding  or  subtracting  the  phase  angle  adjustment  computed 
for  a  given  position  Ax,  from  equation  (A. 25)  to  the 


phase  angle  computed  at  the  original  position  x,  computed 
from  equation  (A. 17); 

<J>  «  <$>  +  (p  .  (A. 26) 

The  value  for  <J>  is  then  substituted  into  equation  (A. 15) 
for<(>#  and  the  wave  reconstructed  with  the  effects  of 
dispersion  accounted  for. 


APPENDIX  B 

THEORY  OF  MEASUREMENT 
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The  appendix  presents  the  derivation  of  the  theory  used  to  reduce 
data  from  split-Hopkinson  bar  experiments.  It  follows  the  development 
of  Kolsky  (1963). 

If  it  is  assumed  that  a  wave  generated  from  a  disturbance  produced 
in  a  isotropic  elastic  bar  propagates  parallel  to  the  longitudinal  axis 
of  the  bar  such  that  it  can  be  considered  to  be  a  plane  wave,  the 

general  form  of  the  wave  equation  will  be  one-dimensional.  In  the  one¬ 
dimensional  case,  the  elements  of  the  bar  will  extend  and  contract  with 
no  lateral  displacement  of  the  axis  of  the  bar  (Kolsky,  1963).  A  one¬ 
dimensional  condition  is  idealized  in  figure  B.l.  The  assumptions  of 

the  one-dimensional  condition  are:  (1)  each  plane  cross  section  will 

remain  plane  while  in  motion,  (2)  the  stress  over  the  cross  section  will 
be  uniform,  and  (3)  that  the  wavelength  of  the  wave  will  be  long 

compared  to  the  cross  sectional  dimension  of  the  bar  (Kolsky,  1963). 

Isolating  a  small  element  of  the  bar,  dx,  with  cross  sectional  area 
A,  the  stress  at  section  A-A  will  be  <j  and  the  stress  at  section  B-B 
will  be  o  +  (3o/3x)  dx.  Taking  the  x-axis  as  the  direction  of  wave 
propagation,  the  components  of  Newton's  second  law  of  motion  (F  =  ma, 
where  F  is  force,  m  is  mass,  and  a  is  acceleration)  can  be  written  as; 


where  p  is  the  mass  density  of  the  bar  and  U  is  displacement.  Combining 
terms  yields; 


* » '  ■  r, 


||  dx  A  -  p  A  dx  0  . 


(B.5) 


Equation  (2.5)  reduces  to 


da  32U 


3x  "  p  ft5 


(B.6) 


Using  Hooke’s  law  for  an  isotropic  elasto  solid  in  uniaxial  stress; 


o  -  Ee  , 


(B.7) 


and  rewriting  e  in  terms  of  displacement  as; 


(B.  8) 


equation  (2.7)  can  be  written  as; 


„  _  E  M 

0  E  3x  • 


(B.9) 


Taking  the  partial  derivative  of  equation  (B.9)  with  respect  to  x 


yields; 


do  „  32U 
3x  ’  3xs 


(B. 10) 


Substituting  the  expression  for  3o/3x  from  equation  (B.10)  into  equation 


(B.6)  yields; 
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(A. 11) 


32U  _  E  32U 
3t2  “  p  3x2  * 


2 

Substituting  Cq  -  E/p  into  equation  (B.11)  provides  an  expression  for 
the  propagation  of  longitudinal  waves  of  infinite  wavelength  along  a 
bar; 


32U  .2  32U 

Jt7  "  0  3X7  * 


(B. 12) 


The  solution  of  equation  (B.12)  is  attributed  to  d'Alembert,  and  may  be 
given  as  (Timoshenko  and  Goodier,  1970); 


U  -  f(x  -  C  t)  +  g(x  +  C  t)  ,  (B. 13) 

o  o 

where  f  and  g  are  arbitrary  functions,  with  f  representing  a  wave 
travelling  in  the  positive  x-direction  and  g  representing  a  wave 
travelling  in  the  negative  x-direction. 

If  a  wave  is  considered  to  travel  only  in  the  negative  x-direction, 
equation  ( B. 13)  c an  be  written  as; 


U  -  g(x  +  CQt)  . 


(B. 1 4) 


Differentiating  equation  (B. 14)  with  respect  t  x  yields; 


2”  .  g, 

3x  8  ’ 


(B.1 5) 


where  the  prime  denotes  the  derivative  with  respect  to  the  argument 
Differentiating  equation  < B .  1  it )  with  respect  to  t  yields; 
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Combining  equations  (B.15)  and  (B.16)  yields; 


3U  3U 

3t  o  3x  * 


(A. 17) 


Substituting  Hooke's  law  (equation  B.9)  into  equation  (B.17)  yields; 


3U  Coq 
3t  "  E 


(B. 1 8) 


Substituting  pC2Q  for  E,  equation  (B.18)  can  be  rewritten  as: 


By  substituting  v  for  the  particle  velocity  (3U/3t),  equation  (B.19) 
simplifies  to; 


a  =  pC  v  .  (B.20) 

o 

The  product  pCQ  is  commonly  referred  to  as  the  characteristic  impedance 
of  the  material  (Rinehart,  1975).  The  procedure  followed  in  deriving 
the  equation  for  particle  velocity  is  similar  to  that  of  other  authors 
(e.g.,  Kolsky,  1963.  Rinehart,  1975;  Zukus,  Nicholas,  Greszczuk,  and 
Curran,  1982). 

A  diagram  of  the  SHPB  near  the  specimen  is  shown  in  figure  B.2. 

The  stress  waves  a.  (incident  wave)  and  o  (relfected  wave)  act  on 

i  r 

interface  1  and  ot  (transmitted  wave)  acts  on  interface  2.  Before  any 


relfections  occur,  the  particle  velocity  of  the  incident  bar  is  given 
as; 


v 


i 


(B.21  ) 


If  the  characteristic  impendance,  or  area  of  the  specimen,  is  different 
than  that  of  the  pressure  bars,  a  portion  of  the  compressive  stress  wave 
at  interface  1  will  be  reflected  as  a  tensile  wave  and  that  portion  of 
the  stress  wave  which  the  specimen  is  able  to  support  is  transmitted 
through  the  specimen.  When  the  portion  of  the  stress  wave  propagating 
through  the  specimen  reaches  interface  2,  the  wave  is  once  again 
partitioned,  with  a  portion  being  transmitted  into  the  transmitter 
bar.  The  reflected  wave  at  interface  2  is  compressive;  hence,  it  will 
continue  to  traverse  the  specimen,  increasing  in  amplitude  with  each 
transit.  The  characteristic  impedance  of  the  specimen,  relative  to  the 
characteristic  impedance  of  the  pressure  bars,  and  any  difference  in 
area  will  govern  the  increase  in  amplitude.  The  result  of  these 
internal  reflections  is  that  the  net  particle  velocity  of  interface  1 
will  increase,  and  then  decrease  with  time,  while  the  net  stress  on 

i 

interface  1  will  decrease,  and  then  increase  with  time.  The  net  j 

particle  velocity  of  interface  2  will  Increase  with  time,  as  will  the  I 

t 

net  stress  on  interface  2.  < 

i 

I 

The  equations  for  the  reflection  and  transmission  coefficients  for  j 

! 

a  wave  impinging  on  an  interface  are  derived  in  figure  B.3  Figure  B.4  i 


illustrates  the  process  of  elastic  wave  propagation  exercising  the 
repeated  application  of  the  relfection  and  transmission  coefficients. 
The  process  is  terminated  at  the  point  of  stress  equilibrium  with  the 


Incident  wave.  this  figure  was  prepared  with  the  assumption  that  the 
wavelength  of  the  incident  wave  is  infinite.  For  illustrative  purposes, 
it  was  also  assumed  that  the  density  of  the  specimen  was  one-half  that 
of  the  pressure  bars  and  the  wave  velocity  ( Cq )  in  the  specimen  and 
pressure  bars  is  equal,  thus  the  characteristic  impendance  of  the 
specimen  is  one-half  that  of  the  pressure  bars.  The  pressure  bars  are 
considered  to  be  of  a  sufficient  length  to  avoid  the  necessity  of 
considering  reflections  from  their  ends.  An  implicit  assumption  for 
using  equations  (4)  and  (5)  derived  in  figure  B.3  is  that  do/de  is 
constant  for  the  specimen  material.  However,  when  a  specimen  deforms 
plastically,  as  do  the  soil  specimens  in  this  research,  do/de  is  not 
constant  but  will  be  dependent  on  the  velocity  at  which  the  wave 
propagates  in  the  specimen. 

The  net  particle  velocity  of  interface  1  will  be; 


vi  -  vi  - V  • 


(B. 22) 


°i  -  (“<Jr) 

pc 


(B. 23) 


and  the  particle  velocity  of  interface  2; 


v 


2 


(B.24) 


Representative  particle  velocities  for  interfaces  1  and  2  with  soil  as 
the  specimen  are  shown  in  figure  B.5.  By  taking  the  difference  of  the 
particle  velocities  at  each  interface,  the  rate  at  which  the  specimen  is 


straining  can  be  computed  as; 


(V  ("  vr)  "  vt 


(B. 25) 


(Of  (-  or)  -  ot) 


(B. 26) 


The  strain  experienced  by  the  specimen  at  any  time  t,  can  be  computed  by 


taking  the  integral  of  the  strain-rate; 


«  J  e  dt  . 
0 


(B.27) 


Figures  B.6  and  B.7  show  typical  results  for  e  and  e;  respectively,  as 


computed  from  equations  (B.26)  and  (B.27).  If  the  specimen  is  straining 


uniformly  then  the  values  obtained  from  equations  (B.26)  and  (B.27)  will 


be  exact.  However,  the  specimen  will  not  usually  experience  uniform 


straining  over  its  length.  Therefore,  the  computed  strain-rate  and 


strain  can  only  be  considered  average  values  for  the  specimen. 


The  stresses  at  interfaces  1  and  2  can  be  determined  in  a  similar 


fashion.  The  force  at  interfaces  1  and  2  will  be; 


(Oi  +  (-  or)) 


°t  A2  • 


(B. 28) 


(°t  +  (~  °r^  k2 


°t  A1  ’ 


(B. 29) 


where  is  the  stress  transmitted  into  the  specimen,  is  the  stress 


reflected  at  interface  2,  A1  is  the  area  of  the  pressure  bars,  and  A2  is 


the  area  of  the  specimen.  The  stress  at  each  interface  is  then; 
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°tA1 

’2  "  a2  * 


(B.30) 


(B.3D 


and  the  average  stress  in  the  specimen  will  be; 


(B.32) 


or; 


(a.  +  (-  ar)  ♦  ot)  A1 
°avg“  2  A2 


(B.33) 


Figure  B.8  shows  a  graphic  representation  of  the  averaging  procedure  and 
its  result. 


9 


tft 


At  interface  A  -  A, 

A 

ir  -  0. 

<0i  4  V  A1  •  °tA2, 

(1) 

Continuity 

vi  '  vr  **  V 

(2) 

From  equation  (2.20), 

equation  (2)  can  be  written  as, 

°i  °r  °t 

P1C1  ‘  Vl  "  P2C2  ' 

(3) 

Solving  equation  (3)  for  ot#  and  substituting  the  result  into  equation  (1), 
a  solution  for  in  terns  of  can  be  expressed  as  follows# 


o 

r 


A  P  C 

in 


Aap2ca  +  Aipici 


o  • 
1 


By  following  the  ease  procedure  end  solving  equation  (3)  for  or,  and 
substituting  the  result  into  equation  (1),  a  solution  for  in  terns  of 
Oj  can  be  expressed  as. 


(4) 


o 

t 


2A1B2C2 


A2P2C2  *  A1P1C1 


15) 


e 


Figure  B.3  Derivation  of  reflection  and  transmission  coefficients 
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APPENDIX  C 

DERIVATION  OF  THE  RELATIONSHIP  BETWEEN 
THE  STRIKER  BAR  IMPACT  VELOCITY  AND 
THE  MAGNITUDE  OF  THE  APPLIED 
STRESS  WAVE 


The  relationship  between  the  striker  bar  impact 
velocity  and  the  magnitude  of  the  applied  stress  wave  is 


accomplished  through  momentum  considerations  derived 
through  Newton's  second  law  of  motion; 

F  =  m  a  ,  (C.l) 

where  F  is  force,  m  is  mass,  and  a  is  acceleration. 

Equation  (C.l)  can  be  rewritten  by  making  the  following 
substitutions ; 


F  =  O  A  , 


(C.2) 


a  *  _dV  , 
dt 


(C.3) 


where  0  is  stress,  A  is  the  cross  sectional  area  of  the 
striker  bar,  V  is  velocity,  and  t  is  time,  as; 


The  limit  of  integration  t,  on  the  left  side  of  equation 
(C.4 )  is  the  duration  of  impact  between  the  striker  bar 
and  the  incident  bar.  When  the  striker  bar  impacts  the 
incident  bar,  a  compressive  wave  is  generated  in  the 
incident  bar  and  striker  bar  as  well.  The  compressive  wave 
in  the  incident  bar  propagates  toward  the  sample.  The 
compressive  wave  in  the  striker  bar  propagates  in  the 
opposite  direction,  toward  the  free  end  of  the  striker 
bar.  Upon  reaching  the  free  end,  the  compressive  wave  is 


reflected  as  a  tensile  wave  which  returns  to  the  impact 
interface.  As  the  interface  is  unable  to  sustain  tension, 
unloading  occurs,  terminating  the  applied  stress  wave. 
Hence,  the  duration  of  impact  is  the  time  required  for  a 
wave  to  traverse  the  length  of  the  striker  bar  twice,  and 
can  be  expressed  as; 


t  *  2  L  ,  (C.5) 

Co 

where  L  is  the  length  of  the  striker  bar,  and  Co  is  the 
bar  velocity  of  the  striker  bar.  As  mentioned  previously, 
a  0.254  in  long  striker  bar  was  used  in  this  research, 
hence  the  duration  of  the  applied  stress  wave  is 
approximately  100  microseconds.  The  limit  of  integration 
V,  on  the  right  side  of  equation  (C.4)  is  the  striker  bar 
impact  velocity.  Carrying  out  the  integration  of  equation 
(C.4)  yields; 

a  A  2  L  *  m  V  .  (C.6) 

Co 

Substituting  p  A  L  *  m  reduces  equation  (C.6)  to; 

O  «PCn  V  .  (C . 7 ) 

2 


The  magnitude  of  the  applied  stress  is  therefore 
determined  by  the  velocity  at  which  the  striker  bar 
impacts  the  incident  bar. 
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EXPERIMENTAL  ASSUMPTIONS  OF  THE  SPLIT 
HOPKINSON  BAR  EXPERIMENT 


-  125  - 


The  limitations  of  the  SHPB  experimental  method  are  dependent  on 
how  well  the  assumptions  required  to  reduce  the  data  are  satisfied. 
These  assumptions  are: 

(1)  there  is  a  uniform  distribution  of  axial  and  radial  stress 
over  the  length  of  the  specimen, 

(2)  the  waves  in  the  pressure  bars  propagate  without  dispersion, 

(3)  the  stress  state  over  the  cross  sectional  area  of  the  pressure 
bars  is  one-dimensional,  and 

(*0  the  interfaces  between  the  pressure  bars  and  the  specimen  are 
frictionless. 

These  are  the  basic  assumptions  made  Kolsky  (19*19).  Each  assumption, 
its  bearing  on  the  experimental  results  and  the  degree  to  which  the 
assumptions  are  satisfied  for  the  experiments  on  soil  specimens  reported 
herein  follows. 

(1)  Uniform  Distribution  of  Stress 

As  the  stress  wave  first  enters  the  specimen  axial  and  radial 
inertia  forces  act  to  oppose  the  equilibration  of  stress.  This 
establishes  a  stress  gradient  between  the  two  faces  of  the  specimen.  If 
the  wave-transit  time  in  the  specimen  is  small  compared  to  the  duration 
of  the  applied  stress  wave,  many  reflections  can  take  place  within  the 
specimen.  The  stress  gradient  will  diminish,  resulting  in  a  nearly 
uniform  distribution  of  stress  over  the  length  of  the  specimen.  If 
there  is  insufficient  time  for  the  stress  in  the  specimen  to 
equilibrate,  the  experiment  is  one  of  wave  propagation  rather  than  one 
of  uniform  stress.  Should  this  be  the  case,  a  constitutive  equation 
would  be  required  to  study  the  phenomena.  This  defeats  the  purpose  of 
the  experiment  as  it  is  the  constitutive  equation  that  is  sought  in  the 
first  place.  Hence,  if  care  is  exercised  in  selecting  specimen  geometry 


and  the  duration  of  the  applied  stress  wave,  the  SHPB  experiment 
provides  a  method  of  obtaining  a  material  relationship  without  the 
disadvantages  of  having  to  undertake  a  wave  propagation  analysis. 


This  assumption  is  perhaps  the  most  difficult  to  satisfy  when  using 
soil  specimens  in  a  SHPB  experiment.  Because  soils  exhibit  nonlinear 
hysteretic  behavior  and  have  low  wave  speeds  (e.g.,  =  300  m/s)  there  may 
be  insufficient  time  for  the  stress  gradient  between  the  two  faces  of 
the  specimen  to  diminish.  Therefore,  to  satisfy  the  uniform  stress 
assumption  the  specimen  aspect  ratio  will  have  to  be  small  and  the 
duration  of  the  applied  stress  wave  long. 

Using  the  criterion  developed  by  Davies  and  Hunter  (1963); 


(D.  1 ) 


an  estimate  can  be  established  as  to  when  stress  uniformity  might  be 
achieved  in  a  SHPB  specimen.  The  inequality  is  conservative  because  its 
derivation  is  based  on  a  lower  permissible  value  for  the  propagation 
velocity  of  waves  in  the  specimen  (Davies  and  Hunter, 1963).  Because 
unloading  waves  generated  in  the  experiment  will  travel  faster  than  the 
plastic  waves  in  the  specimen,  quasi-equilibrium  will  be  reached  sooner 
than  indicated  by  the  inequality  (Davies  and  Hunter,  1963). 

Figures  D.1,  D.2,  and  D.3  show  results  of  experiment  1 31* •  Using 
equation  (D.1),  the  stress-strain  curve  (figure  D.1),  the  appropriate 
specimen  length  (1.259  cm),  density  (2.12  g/cm^),  and  T  (=130 

microseconds)  stress  equilibrium  is  estimated  to  be  reached  at  42 
microseconds.  The  values  for  specimen  length  and  density  are  given  in 
table  E.1  of  appendix  E  and  T  is  approximated  from  figure  4.2. 
Numerical  values  for  the  slope  of  the  stress-strain  curve  for  experiment 


13*»  are  given  in  the  example  computer  output  in  appendix  G.  Figure  D.2 
shows  the  stress  histories  at  the  specimen- oar  interfaces.  From  this 
figure  it  can  be  seen  that  after  about  25  microseconds  the  stresses  at 
the  two  interfaces  begin  to  converge  to  approximately  the  same  value  and 
are  indeed  equal  at  50  microseconds.  The  stress  difference  across  the 
specimen  as  a  function  of  time  and  the  average  stress-time  response  is 
shown  in  figure  D.3.  If  inertia  forces  were  absent,  the  stress 
difference,  would  be  zero.  It  can  be  seen  that  after  25  microseconds 
the  stress  difference  is  very  small  relative  to  the  average  stress  from 
which  it  can  be  concluded  that  the  stress  gradient  in  the  specimen  has 
diminished. 

The  experimental  evidence  presented  in  figures  D.2  and  D.3  support 
the  estimate  computed  from  the  Davies  and  Hunter  (1963)  criterion  as  to 
when  stress  uniformity  is  achieved  in  a  SHPB  soil  specimen.  Hence,  the 
Davies  and  Hunter  criterion  can  be  used  to  estimate  when  stress 
uniformity  is  likely  to  be  achieved  in  a  soil  specimen. 

In  contrast  to  the  criterion  of  Davies  and  Hunter  (1963)  which 
provided  a  reasonalbe  estimate  to  the  time  required  for  the  specimen  to 
reach  equilibrium,  the  criterion  of  Bertholf  and  Karnes  (197*0  is 
violated  in  every  experiment  by  at  least  a  factor  of  2  (see  section 
2.5).  For  example,  experiment  1 3 had  a  maximum  strain  rate  of  800 
s  ^ .  To  satisfy  the  Bertholf  and  Karnes  (197*0  criterion  the  specimen 
strain-rate  and  loading  wave  rise-time  would  have  to  be  less  than 
800  S~1  and  greater  than  90  microseconds;  respectively. 

The  purpose  of  the  Bertholf  and  Karnes  (197*0  criterion  was  to 
bound  the  experimental  parameters  such  that  the  influence  of  inertia 
would  be  minimized.  Although  the  experiments  in  this  research  did  not 


satisfy  the  criterion,  it  has  been  shown  that  stress  equilibrium  will  be 
achieved,  satisfying  the  uniform  stress  assumption. 


The  influence  of  inertia  and  its  affect  on  the  ability  of  a  soil  j 

specimen  to  reach  stress  equilibrium  has  been  investigated  j 

experimentally,  and  the  results  compared  with  the  estimates  computed  j 

from  the  criterion  of  Davies  and  Hunter  (1963)  and  Bertholf  and  Karnes  i 

(1974).  Experimentally  it  has  been  demonstrated  that  the  stress  at  the  > 

i 

specimen-bar  interfaces  will  be  approximately  equal  within  the  rise-time 
of  the  incident  stress  wave.  This  means  that  the  first  1  to  2  percent 
strain  the  specimen  sustains  must  be  disregarded  as  stress  equilibrium 
within  the  specimen  has  not  been  achieved.  It  has  been  shown  that  the 

J 

criterion  of  Davies  and  Hunter  provides  a  reasonable  estimate  as  to  when  j 

stress  uniformity  is  achieved  in  the  specimen  by  direct  comparison  to  j 

i 

the  experimental  results.  In  contrast  to  the  criterion  of  Davies  and  j 

Hunter,  the  bounds  on  the  experimental  parameters  established  by  \ 

Bertholf  and  Karnes  are  too  restrictive  when  applied  to  soil  specimens. 

(2)  Wave  Dispersion 

The  stress  wave  initiated  by  the  impact  of  the  striker  bar  on  the 
incident  pressure  bar  is  initially  very  complex.  This  is  due  in  part  to 
friction  at  the  striker  bar-incident  bar  interface,  as  well  as  the 
appearance  of  spherical  and  dilatational  waves  (Follansbee  and  Frantz, 

1983).  These  end  effects  disappear  at  about  10  bar  diameters  from  the 
point  of  impact  (Wasley,  1973.  and  Yeung  Wye  Kong  et  al.,  1974). 

Figure  D.4  shows  the  waves  as  recorded  by  the  incident  bar  strain 


gauge  in  a  test  in  which  no  specimen  was  placed  between  the  incident  and 
transmitter  bars.  The  incident  and  transmitter  bars  were  separated  such 
that  the  specimen-incident  bar  interface  was  a  free  end.  In  this 


configuration,  the  applied  stress  wave  upon  reaching  the  free  end  is 
reflected  as  a  tensile  wave  with  amplitude  equal,  but  opposite  in  sign 
to  the  applied  wave.  If  the  wave  does  not  disperse  during  propagation, 
the  two  waves  should  be  identical.  Clearly,  wave  dispersion  does  occur 
as  the  period  of  the  oscillations  of  the  reflected  wave  are  greater  than 
those  of  the  incident  wave. 

Figures  D.5  and  D.6  demonstrate  the  application  of  the  dispersion 
correction  (see  appendix  A)  to  the  waves  generated  in  the  SHPB  test  with 
the  specimen-incident  bar  interface  as  a  free  end.  Figure  D.5  compares 
the  reflected  stress  wave  as  recorded  at  the  strain  gauge  location  with 
the  dispersion  corrected  wave  as  reconstructed  at  the  specimen-incident 
bar  interface.  It  is  evident  that  the  dispersion  correction  has 
effectively  compressed  the  reflected  wave  during  the  reconstruction. 
This  is  to  be  expected  as  the  phase  angle  will  decrease  for  negative 
values  of  Ax  (see  appendix  A,  equation  A. 26).  In  figure  D.6  the 
dispersion  corrected  incident  and  reflected  stress  waves  as 
reconstructed  at  the  specimen-incident  bar  interface  are  shown.  Figure 
D.6  shows  there  is  a  correspondence  of  the  peaks  and  valleys  of  the  two 
waves.  This  correspondence  was  not  evident  between  the  un corrected 
stress  waves  at  the  recorded  strain  gauge  position  (see  figure  D.4). 
This  correspondence  supports  the  assumption  that  only  the  fundamental 
mode  of  vibration  is  excited  in  the  SHPB  experiment  (Follansbee  and 
Frantz,  1983).  If  other  modes  of  vibration  had  been  present,  the 
procedure  of  performing  the  dispersion  correction  based  on  obtaining  the 
waves'  component  phase  velocity  and  wavelength  from  the  fundamental  mode 
would  not  have  been  successful  (Follansbee  and  Frantz,  1983). 

Figure  D.iJ  clearly  shows  that  wave  dispersion  does  occur,  while 


figure  D.6  shows  that  the  Follansbee  and  Frantz  (1983)  correction 
procedure  can  be  used  to  account  for  the  phenomenon.  If  dispersion  is 
not  accounted  for  in  the  reduced  experimental  data,  large  oscillations 
will  appear  in  the  stress-strain  curve  which  can  mask  trends  and 
increase  difficulty  in  interpreting  specimen  response. 

Davies  (1948)  showed  that  oscillations  in  the  recorded  stress  waves 
were  attributable  to  dispersion  and  could  be  predicted.  However,  until 
recently  there  has  been  no  technique  that  could  easily  be  incorporated 
into  the  standard  SHPB  data  reduction  procedure  to  account  for  this 
phenomenon  (Follansbee  and  Frantz,  1983).  Consequently,  the 
oscillations  in  the  raw  data  were  retained  and  a  smoothing  technique 
applied  to  the  reduced  data.  Although  the  smoothing  of  reduced  SHPB 
data  is  an  accepted  practice  (e.g.,  Rand,  1967),  a  systematic  method  of 
accounting  for  the  dispersive  nature  of  a  propagating  wave  is 
preferable. 

The  effects  of  dispersion  on  the  experimental  results  obtained  in 
this  research  was  investigated  by  conducting  a  test  with  the  end  of  the 
incident  bar  being  a  free  end.  The  correction  procedure  was  shown  to 
account  for  the  wave  dispersion  occuring  between  the  strain  gauge 
location  and  the  free  end  of  the  incident  bar. 

Figure  D.7  illustrates  the  effect  of  dispersion  on  the  appearance 
of  the  stress-strain  curve  for  experiment  145.  Here  it  can  be  seen  that 
the  dispersion  corrected  curve  is  smoother  and  easier  to  interpret.  For 
instance,  it  would  have  been  difficult  to  correlate  the  sudden  increase 
in  stiffness  to  the  condition  of  zero  air  voids  had  the  dispersion 
correction  procedure  no  been  employed.  Hence,  an  effect  of  dispersion 
on  experimental  results  is  to  mask  trends. 


Stress  Variation  Over  the  Cross  Section  of  the  Bar 


The  information  required  to  evaluate  the  strain-rate-time,  strain¬ 
time,  and  stress- time  response  of  a  specimen  in  a  SHPB  experiment  is 
obtained  from  strain  gauges  mounted  on  the  surface  of  the  pressure 
bars.  If  this  recorded  information  is  not  consistent  with  the  stress¬ 
time  response  the  pressure  bar  is  experiencing  over  its  cross  sectional 
area,  the  resulting  computations  may  be  seriously  in  error.  Hence,  a 
significant  condition  that  must  be  present  in  any  SHPB  experiment  is 
that  the  longitudinal  stress  and  displacement  be  nearly  constant  over 
the  cross  sectional  area  of  the  bar. 

A  prediction  of  the  dispersion  equation  is  that  the  longitudinal 
stress  and  displacement  will  vary  over  the  cross  sectional  area  of  the 
pressure  bars,  and  that  the  variation  will  be  dependent  on  the  ratio 
R/A  (Davies,  19*18).  Therefore,  errors  are  unavoidable  if  one¬ 

dimensional  wave  analysis  is  used  to  reduce  the  data.  The  significance 
of  the  error  and  its  impact  on  interpreting  experimental  results  was 
investigated  by  Davies  (19*18).  Using  the  limiting  condition  of  R/A  «  0 
as  an  error  baseline  (i.e.,  theoretically  at  R/A  -  0  the  longitudinal 
displacement  and  stress  are  constant  over  the  cross  sectional  area  of 
the  bar),  Davies  was  able  to  establish  that  for  R/A  <  0.1,  the 
displacement  recorded  at  the  surface  of  the  pressure  bar  would  differ 
from  that  at  the  bars  longitudinal  axis  by  ±  3  percent  (Davies,  19*i8). 

In  this  research  the  diameter  of  the  pressure  bars  are  larger  than 
those  used  by  past  investigators  by  a  least  a  factor  of  2  and  In  some 
cases  a  factor  of  3.  For  example,  in  the  experiments  performed  by 
Lindholm  (196*1)  the  diameter  of  the  pressure  bars  were  1.9  cm.  A 


concern  in  using  the  larger  bars  is  that  the  behavior  over  the  cross 
section  will  deviate  from  one-dimensional  condition  that  is  assumed; 
hence  increasing  the  error  in  the  computations  beyond  that  established 
by  Davies. 

The  diameter  of  the  pressure  bars  used  in  this  research  was 
significantly  larger  than  that  used  in  past  SHPB  experimentation.  A 
consequence  of  using  a  larger  diameter  bar  is  that  the  stress-time 
history  over  the  cross  sectional  area  of  the  pressure  bar  may  not  be 
constant.  Hence,  measurements  made  at  the  surface  of  the  bar  might  not 
be  representative  of  the  behavior  over  its  cross  section.  In  the  work 
of  Davies  (1948)  he  showed  that  if  R/A  <  0.1  the  displacements  at  the 
bar  surface  would  differ  from  those  along  the  axis  of  the  bar  by  ±  3 
percent . 

In  order  to  investigate  the  degree  to  which  the  criterion 
established  by  Davies  is  satisfied  for  the  SHPB  used  in  this  research, 
the  incident  wave  for  experiment  134  was  decomposed  into  its  respective 
Fourier  components.  Table  D.1  lists  the  values  n,  R/A,  Cn/CQ,  and  the 
ratios  of  the  amplitudes  of  the  Fourier  component  nu>o(Dn)  to  the  largest 
magnitude  Fourier  component  (AQ). 

For  n  <  9,  the  value  R/A  is  less  than  0.1,  and  the  magnitude  of  the 
Fourier  component  (DR)  is  approximately  6  percent  of  the  largest 
magnitude  Fourier  component  (AQ).  These  computations  illustrate  that 
the  high  frequency  Fourier  components  of  the  wave  are  damped  out  quickly 
which  tends  to  support  a  nearly  uniform  stress-time  distribution  over 
the  cross  sectional  area  of  the  bar.  However,  n  -  17  is  required  to 
model  the  experimentally  observed  wave;  therefore,  the  error  in  the 
computations  will  be  in  excess  of  ±  3  percent. 


Figure  D.8  shows  the  incident  wave  modeled  with  n  -  17,  and  n  -  9 
along  with  their  difference.  The  difference  was  taken  as  a  qualitative 
indicator  of  the  error  in  excess  of  ±  3  percent.  A  quantitative  measure 
can  be  computed  by  summing  the  absolute  values  of  the  Fourier 
coefficients  for  n  equal  10  to  17  and  dividing  by  the  summation  of  the 
Fourier  coefficients  for  n  equal  0  to  17. 
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This  relative  error  will  be  an  upper  bound  as  Dn  is  computed  from  An  and 
Bn  which  can  be  opposite  in  sign.  Hence,  as  a  worst  case  the 
computational  error  will  be  approximately  ±  10  percent. 

Figure  D.9  shows  the  difference  in  the  computed  average  stress- 
strain  response  between  using  n  -  17  and  n  -  9.  It  can  be  seen  that  the 
stress-strain  response  is  not  significantly  altered;  hence  ±  10  percent 
can  be  considered  to  be  an  upper  bound  on  the  error  when  assuming  one¬ 
dimensional  wave  analysis  to  reduce  the  data  obtained  using  60.3  mm 
diameter  pressure  bars. 


Specimen-Bar  Interface  Friction 

In  the  traditional  SHPB  experiment,  a  specimen  with  a  diameter 
slightly  less  than  that  of  the  bars  is  placed  between  them  (Lindholm, 
1 964 ) .  This  is  to  allow  the  specimen  to  expand  radially  during  the 
experient  while  not  exceeding  the  diameter  of  the  pressure  bars.  As 
stress  is  applied  to  the  specimen,  radial  shear  stresses  are  created 
between  the  pressure  bars  and  the  specimen.  This  has  been  commonly 


referred  to  as  "end  effects"  or  "friction  effects."  A  result  of  these 
effects  is  that  the  specimen  tends  to  be  clamped  at  the  specimen-bar 
interfaces  (Green  and  Perkins,  1969).  This  prevents  the  specimen  from 
expanding  uniformly;  hence,  barreling  of  the  specimen  is  often 
observed.  This  may  result  in  a  strengthening  effect  being  observed  in 
the  stress-strain  response  may  greatly  increase  the  degree  of 
nonuniformity  of  stress  and  strain  in  the  specimen  (Bertholf  and  Karnes, 
197*0. 

To  reduce  the  effects  of  friction  it  is  common  practice  to  apply  a 
film  of  lubrication  to  the  ends  of  the  pressure  bars  in  contact  with  the 
specimen.  The  results  of  the  numerical  study  conducted  by  Bertholf  and 
Karnes  (197*1)  have  shown  this  to  be  an  acceptable  method  to  minimize  the 
effects  of  friction  on  stress-strain  response.  However,  the  effects  of 
friction  cannot  be  considered  independently  of  specimen  geometry,  i.e., 
aspect  ratio  (Rand,  1967).  Rand  (1967)  has  shown  that  decreasing  the 
specimen  aspect  ratio  will  have  the  same  effect  as  increasing  the 
coefficient  of  friction  between  the  specimen  and  the  pressure  bars. 
Again  as  a  result  of  the  workof  Bertholf  and  Karnes  (197*0,  a  specimen 
aspect  ratio  of  0.5  is  thought  to  be  an  optimum. 

However,  for  soil  specimens  an  aspect  ratio  of  0.5  may  be  too 
restrictive  due  to  the  same  concerns  mentioned  in  section  D.1  (i.e.,  low 
wave  speed  and  nonlinear  hysteretic  behavior).  Therefore  the  aspect 
ratio  f  >r  a  soil  specimen  msut  be  large  enough  so  that  friction  effects 
are  minimized  but  small  enough  to  permit  stress  equilibration  within  the 
duration  of  the  applied  stress  wave.  A  method  to  retain  a  small  aspect 
ratio  while  avoiding  friction  effects  is  to  contain  a  soil  specimen  in  a 
confining  cylinder.  The  confining  cylinder  prevents  barreling  and 


excessive  radial  strain  so  that  the  specimen  is  in  a  state  of  nearly 
uniaxial  strain.  Hence,  friction  effects  on  specimen  stress-strain 
response  is  minimized. 

The  numerical  work  of  Bertholf  and  Karnes  (1974)  has  shown  that 
friction  between  the  specimen  and  the  pressure  bars  can  affect  specimen 
response,  particularly  for  relatively  thin  specimens.  Interface 
friction  affects  specimens  response  in  two  ways:  (1)  a  greater  stress 
will  be  required  to  produce  a  given  strain,  and  (2)  the  specimen  will 
barrel.  From  their  work  they  concluded  that  the  effects  of  friction 
could  be  minimized  if  the  specimen-bar  interfaces  were  lubricated  and  a 
specimen  aspect  ratio  of  approximately  0.5  selected. 

All  the  specimens  in  this  research  were  the  same  diameter  as  the 
pressure  bars  and  placed  in  a  confining  cylinder  which  prevented 
excessive  barreling  and  hence  negligible  shear  stresses  between  the 
specimen  and  the  pressure  bars.  This  was  verified  in  section  5.3  where 
it  was  shown  that  the  specimen  did  sustain  some  radial  strain  but  for 
all  experiments  was  less  than  1  percent. 

The  results  of  chapter  5  can  be  used  to  examine  the  effect  of 
friction  on  specimen  stress-strain  response.  In  figures  5.11,  5.12,  and 
5.13  the  stress-strain  response  is  compared  for  different  length 
specimens  to  the  nominal  applied  stresses  of  250,  400,  and  520  MPa; 
respectively.  These  figures  are  useful  because  the  specimen  lengths  are 
different  by  a  factor  of  two;  hence  if  friction  effects  exist  their 
stress-strain  response  should  not  be  the  same. 

In  all  three  figures  the  stress-strain  response  for  the  different 
length  specimens  is  nearly  indent! cal  up  to  their  initial  gas 
porosities.  For  strains  in  excess  of  the  initial  gas  porosity  the 


specimen  response  is  governed  principally  by  the  pore  water.  This  is 
evident  from  the  results  presented  throughout  chapter  5.  In  all  three 


r 


. 


figures  the  peak  stress  of  the  shorter  specimens  exceeded  that  of  the 
longer  by  at  least  25  percent.  As  the  specimen  does  sustain  some  radial 
strain  this  phenomenon  may  be  due  to  friction.  However,  to  be 
conclusive,  experiments  need  to  be  conducted  on  specimens  with  as 
similar  initial  gas  porosities  as  possible. 

Hence,  for  soil  specimens  with  lengths  different  by  a  factor  of  two 
and  aspect  ratios  less  than  0.2,  friction  is  not  observed  to  influence 
the  stress-strain  response  at  strains  less  than  the  initial  gas 
porosity.  From  this  conclusion  it  can  be  determined  that  the  optimum 
aspect  ratio  of  0.5  is  too  restrictive  when  applied  to  soil  specimens 
and  that  by  maintaining  a  condition  of  nearly  uniaxial  strain,  aspect 
ratios  less  than  0.5  can  be  used  with  the  effects  frction  on  stress- 
strain  response  held  to  a  minimum. 


R/A 


.Cn/Co 


Dn/^o 


0.01170 

0.02340 

0.03500 

0.04670 

0.05840 

0.07020 

0.08200 

0.09380 

0.10570 

0.11770 

0.12980 

0.14210 

0.15450 

0.16730 

0.18030 

0.19380 

0.20770 

0.22220 

0.23750 

0.25370 


0.99966 

0.99918 

0.99864 

0.99804 

0.99733 

0.99646 

0.99538 

0.99402 

0.99229 

0.99012 

0.98743 

0.98408 

0.98002 

0.97505 

0.96911 

0.96196 

0.95353 

0.94358 

0.93186 

0.91817 


0.89275 

0.65881 

0.41980 

0.18439 

0.02868 

0.12833 

0.13652 

0.11074 

0.05668 

0.02368 

0.05596 

0.05409 

0.04164 

0.02027 

0.01707 

0.02908 

0.02016 

0.00962 

0.01329 

0.00912 
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SUPPORTING  TABLES 
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Super  iaent 

Specimen 

Moisture 

Gas 

Met 

Void 

Degree  of  j 

Mo. 

length 

Content 

Porosity 

Density 

Ratio 

Saturation 

(OB) 

(%> 

(%) 

(g/oc  ) 

1 

«** v 


i 


A 


Table  E.l  cont. 


Experiment 

MO. 

Specimen 

Length 

(o» 

Moisture 

Content 

(«) 

Gee 

Porosity 

(%> 

Net 

Density 
(g/oc  ) 

Void 

Ketio 

Degree  of 
Setureticn 
(») 

146 

1.267 

11.9 

6.26 

2.13 

0.41 

78.3 

147 

0.636 

11.9 

6.64 

2.12 

0.41 

77.2 

148 

0.631 

14.4 

7.41 

2.04 

0.50 

77.6 

162 

0.629 

13.0 

4.44 

2.14 

0.41 

•4.7 

163 

0.631 

12.9 

4.84 

2.13 

0.41 

83.4 

164 

1.298 

12.9 

7.49 

2.07 

0.45 

76.0 

165 

1.289 

12.7 

7.05 

2.09 

0.44 

77.0 

166 

0.624 

14.0 

4.42 

2.12 

0.44 

85.5 

167 

0.622 

14.0 

4.15 

2.12 

0.43 

•6.3 

Table  E.2 


Moisture  Loss  Between  the  Time  a  Specimen  is  Prepared 
and  the  Time  it  is  Positioned  Between  the  Bars 


Experiment  No. 

MCI (a) 

(%) 

MC2 ( b) 
(%) 

Percent  Difference 
(-> 

A 

13.4 

12.3 

8.2 

B 

13.0 

12.2 

6.2 

C 

12.9 

12.2 

5.4 

D 

13.0 

12.7 

2.3 

E 

13.1 

12.6 

3.8 

P 

13.2 

12.6 

4.5 

G 

13.2 

12.4 

6.1 

H 

13.3 

12.4 

6.8 

Arithmetic  mean  *  5.4 

Standard  Deviation  *  1.8 

a)  MCI  is  the  specimen  moisture  content  at  the  time  it  was 
prepared. 

b)  MC2  is  the  specimen  moisture  content  at  the  time  it  was 
positioned  between  the  pressure  bars. 

c)  The  specimens  for  these  tests  were  prepared  following 
the  procedure  that  was  outlined  in  section  3.3.  The 
same  sequence  of  steps  from  specimen  preparation  to 
just  before  firing  of  the  gas  gun  were  followed  as  if 
an  actual  experiment  were  to  be  performed. 
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Table  E.3 


Experimental  Results 
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Table  E.3  cont 


Experiment  No. 

Applied 

Stress 

(MPa) 

Peak.  Stress 

(MPa) 

Strain  at 

Peak  Stress 

(«) 

Strairrtlate  at 

Peak  Stress 

(  s1) 

148 

249 

117 

15.2 

972 

162 

519 

325 

20.6 

1327 

163 

522 

425 

15.3 

459 

164 

507 

309 

12.4 

696 

165 

254 

76 

9.3 

136 

166 

261 

124 

13.3 

967 

167 

260 

151 

12.6 

779 

a)  In  experiment  63  the  gain  netting  for  the  transmitter  bar  strain  gauge  wi 
incorrectly  set.  This  caused  the  digitizer*  collecting  the  strain  gauge  data  to  be 
overdriven.  As  a  result,  only  two  of  the  required  three  data  sets  for  data  reduction 
eere  collected. 

b)  In  experiment  70  a  tight  fit  between  all  system  ooofonents  was  not  Maintained.  As  a 
results,  an  irregular  incident  stress  wave  tats  initiated  that  Mas  wisuitable  for  data 
reduction. 
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Table  E. 4 


Specimen  Seating  Strain 


Experiment 

No. 

Preplacement 

Length 

(cm) 

In- Position 
Length 
(cm) 

Seating 

Strain 

(t) 

112 

0.673 

0.644 

4.3 

113 

0.674 

0.645 

4.3 

114 

0.723 

0.655 

9.4 

115 

0.677 

0.645 

4.7 

116 

0.715 

0.637 

10.9 

117 

0.697 

0.635 

8.9 

118 

0.694 

0.645 

7.1 

119 

0.721 

0.645 

10.5 

131 

1.349 

1.313 

2.7 

132 

1.338 

1.269 

5.2 

133 

1.367 

1.291 

5.6 

134 

1.382 

1.259 

8.9 

135 

1.337 

1.307 

2.2 

136 

1.330 

1.223 

8.0 

137 

0.670 

0.665 

0.7 

138 

0.688 

0.596 

13.4 

139 

0.728 

0.609 

16.3 

145 

1.360 

1.265 

7.0 

146 

1.316 

1.267 

3.7 

147 

0.759 

0.636 

16.2 

162 

0.646 

0.629 

2.6 

163 

0.684 

0.631 

7.7 

164 

1.362 

1.298 

4.7 

165 

1.345 

1.289 

4.2 

166 

0.687 

0.624 

9.2 

167 

0.654 

0.622 

4.9 

A 

0.727 

0.641 

11.8 

B 

0.711 

0.623 

12.4 

C 

0.715 

0.631 

11.7 

D 

1.395 

1.273 

8.7 

E 

1.364 

1.275 

6.5 

F 

1.369 

1.273 

7.0 

G 

1.335 

1.273 

4.6 

Arithmetic  mean  *  7.5 
Standard  Deviation  *  3.9 

a)  Experiments  numbers  A  through  G  do  not  represent 
complete  experiments.  These  were  tests  performed  to 
provide  additional  information  on  seating  strains  and 
moisture  loss. 
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•  Table  E.5 

Computed  Radial  Displacements  Due  to  Applied  Load 


1 

j# 

1 

i 

Experiment  No. 

Inside  Radius  of  Cylinder 
( cm) 

Displacement 

(cm) 

> 

j 

31 

3.017 

0.0007 

32 

3.018 

0.0264 

38 

3.025 

0.0191 

39 

3.030 

0.0094 

40 

3.059 

0.0120 

i 

41 

3.055 

0.0147 

42 

3.018 

0.0037 

55 

3.048 

0.0177 

• 

56 

3.030 

0.0154 

57 

3.058 

0.0093 

58 

3.055 

0.0064 

59 

3.058 

0.0251 

60 

3.028 

0.0261 

61 

3.033 

0.0118 

• 

62 

3.020 

0.0172 

63 

3.045 

****** 

70 

3.030 

****** 

71 

3.032 

0.0240 

72 

3.063 

0.0118 

73 

3.035 

0.0201 

• 

74 

3.053 

0.0238 

75 

3.038 

0.0082 

76 

3.018 

0.0020 

112 

3.060 

0.0182 

113 

3.056 

0.0177  ! 

114 

3.026 

0.0175 

• 

115 

3.060 

0.0187 

116 

3.060 

0.0111 

117 

3.056 

0.0149 

118 

3.059 

0.0115 

119 

3.040 

0.0109 

131 

3.045 

0.0138  $ 

V 

132 

3.062 

0.0146 

133 

3.060 

0.0124 

134 

3.055 

0.0150  > 

135 

3.030 

0.0025  i 

136 

3.020 

0.0198  * 

137 

3.018 

0.0166  r 

m 

138 

3.025 

0.0129  j 

139 

3.025 

0.0080  * 

145 

3.025 

0.0073  3 

146 

3.028 

0.0220  i 

-  155 


Table  E.5  cont. 


Experiment  Mo. 

Inside  Radius  of  Cylinder 
(cm) 

Displacement 

(cm) 

147 

3.056 

0.0105 

148 

3.035 

0.0083 

162 

3.028 

0.0228 

163 

3.023 

0.0297 

164 

3.025 

0.0217 

165 

3.055 

0.0053 

166 

3.035 

0.0087 

167 

3.028 

0.0106 

Table  E.6 


Mass  of  Soil  Lost  During  the  Experiment 


Experiment  No. 

Mass  Lost 

<g> 

138 

1.7 

145 

2.1 

147 

1.9 

148 

5.3 

162 

7.9 

164 

4.8 

165 

2.1 

166 

4.8 

167 

2.3 

Table  E.7 


Measured  Moisture  Content  Changes  Before 
and  After  the  Experiment 


Experiment  No. 


Preexperiment 
Moisture  Content 


Postexperiment 
Moisture  Content 


Percent 

Difference 


• 

(%) 

(%) 

(- 

31 

11.9 

10.0 

16.0 

32 

11.9 

10.8 

9.2 

38 

13.2 

12.2 

7.6 

• 

39 

13.2 

12.9 

2.3 

40 

13.5 

13.1 

3.0 

41 

13.3 

10.7 

19.5 

42 

13.4 

9.7 

27.6 

55 

12.4 

**** 

**** 

56 

12.6 

**** 

**** 

• 

57 

12.7 

**** 

**** 

58 

13.9 

**** 

**** 

59 

12.7 

**** 

**** 

60 

12.8 

**** 

**** 

61 

13.1 

**** 

**** 

62 

12.9 

**** 

**** 

4» 

63 

13.3 

**** 

**** 

70 

12.1 

10.3 

14.9 

71 

11.7 

9.2 

21.4 

72 

15.4 

8.5 

44.8 

73 

12.2 

10.5 

13.9 

74 

12.0 

11.8 

1.7 

# 

75 

12.3 

10.0 

18.7 

76 

11.9 

10.3 

13.4 

112 

11.8 

10.4 

11.9 

113 

12.1 

10.0 

17.4 

114 

12.1 

10.1 

16.5 

115 

11.4 

10.2 

10.5 

• 

116 

11.2 

9.8 

12.5 

117 

10.7 

9.9 

7.5 

118 

10.6 

9.6 

9.4 

119 

10.4 

9.5 

8.7 

131 

12.4 

11.4 

8.1 

132 

12.5 

11.7 

6.4 

133 

12.5 

11.0 

12.0 

134 

12.4 

11.4 

8.1 

135 

7.0 

6.2 

11.4 

136 

15.1 

11.5 

23.8 

137 

15.4 

12.1 

21.4 

138 

13.0 

12.0 

7.7 

• 

139 

7.0 

6.3 

10.0 

145 

11.8 

11. 3(  4.6) 

4.2 

146 

11.9 

11. 1<  6.2) 

6.7 

147 

11.9 

10. 7(  6.4) 

10.1 
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Table  E.7  cont 


E^ieriment  No. 

Preexperiment 

Postexperiment 

Percent 

Moisture  Content 

Moisture  Content 

Difference 

(«) 

(1) 

(-) 

148 

14.4 

12.8(11.4) 

11.1 

162 

13.0 

11. 7(  6.6) 

10.0 

163 

12.9 

11.9 

7.8 

S- 

164 

12.9 

11. 9(  7.4) 

7.8 

165 

12.7 

12.1(11.1) 

4.7 

■/' 

166 

14.0 

12.7(11.5) 

9.3 

167 

14.0 

12.7(11.0) 

9.3 

v  — 

a)  The  numbers  in  parentheses  beside  the  post  experiment 
moisture  content  percentages  are  moisture  contents  (in 
percent)  of  the  material  scraped  from  the  region  where  the 
bars  and  confining  cylinder  overlap  (see  figure  3.8). 


APPENDIX  P 

PROCEDURE  AND  RESULTS  OF  MOISTURE/DENSITY 
VARIATION  STUDY 


The  soil  used  for  this  study  was  the  same  as  that 
used  in  the  SHPB  experiments.  The  specimens  were  prepared 
at  moisture/density  conditions  near  optimum  as  determined 
from  the  Harvard  miniature  compaction  procedure  (see 
figure  3.5).  The  specimens  for  this  study  were  prepared 
following  the  same  procedures  outlined  in  section  3.3. 

The  specimens  were  compacted  in  stainless  steel 
cylinders  with  the  nominal  dimensions  of  12  cm  in  length 
and  6  cm  inside  diameter.  Static  compaction  of  the 
specimen  was  performed  from  both  ends  of  the  cylinder. 

Four  lifts  of  soil  were  used  to  complete  the  compaction 
processs.  After  compaction,  the  specimen  was  extruded  from 
the  cylinder  and  sliced  into  1  cm  sections. 

The  density  variation  over  the  specimen  length  was 
determined  by  immersing  a  piece  of  each  1  cm  section  in  a 
known  volume  of  mercury,  and  measuring  the  volume  of 
mercury  the  specimen  displaced.  Moisture  content  was 
determined  according  to  the  standard  method  (ASW 
D-2216-80).  The  results  of  the  study  are  presented  in 
figures  F.l  through  F.6. 
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APPENDIX  G 
REDUCE  CODE 
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mm 


tka  winrt  ma  v  a. 


cv  Xv  v 


G.l  Flowchart  for  REDUCE 


reconstruct  wave 
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G . 2  Program  Listing 


REDUCE  CODE  LISTING 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

for 

c 

c 

c 

c 

c 

c 

c*** 

c* 

c* 

c* 

c*** 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


This  program  is  used  to  process  the  data  obtained 
from  a  split-Hopkinson  pressure  bar  experiment.  Two 
sources  of  input  data  are  required. 

1)  The  incident  and  transmitter  bar  strain  gauge 
records  hbxxxi  and  hbxxxt  and  the  calibration 
records  calxxi  and  calxxt.  These  files  will  be 
read  in  from  storage. 

2)  The  experiment  number,  gain  settings,  incident 
wave  start  time,  initial  specimen  length,  and  the 
number  of  Fourier  terms  to  be  used  in  the 
computations.  These  data  are  supplied  manually  by 
the  operator. 

The  output  of  the  program  is  the  average  specimen 
strain-rate-time,  strain-time,  and  stress-time 
histories.  In  addition,  the  particle  velocities  at 
specimen-bar  interfaces  and  stress  difference  between 
the  interfaces  are  computed.  These  data  are  formated 

output  to  a  printer.  The  results  are  also  plotted 
using  the  DISSPLA  9.0  graphics  package.  The 
graphic  subroutines  and  their  calls  have  been  removed 
from  this  listing. 


******************************************************* 

* 

MAIN  PROGRAM  PARAMETER  IDENTIFICATION  * 

* 

******************************************************* 


ARRAYS 


cali  - 

calt  - 

hbi  - 
hbt  - 
barinc  - 


incident  bar  strain  gauge  calibration  storage 
file 

transmitter  bar  strain  gauge  calibration 
storage  file 

incident  bar  strain  gauge  data  storage  file 
transmitter  bar  strain  gauge  data  storage  file 
incident  and  reflected  wave  stress-time 
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c 

trsbar 

- 

c 

diff 

- 

c 

sum 

- 

c 

xint 

- 

c 

xrate 

- 

c 

xstres 

- 

c 

barref 

- 

c 

velin 

- 

c 

c 

velout 

- 

c 

c 

stsdel 

- 

c 

c 

c 

c 

hopout 

- 

c 

hopxxi 

- 

c 

hopxxf 

- 

c 

hopxxt 

- 

c 

c 

hopxxif 

- 

c 

hopxxrf 

- 

c 

hopxxt f 

- 

c 

straxx 

- 

c 

strtxx 

- 

c 

strsxx 

- 

c 

sssrxx 

- 

c 

c 

c 

c 

gaini 

- 

c 

gaint 

- 

c 

sttime 

- 

c 

srtime 

- 

c 

tstime 

- 

c 

test 

- 

c 

expnam 

- 

c 

scfaci 

- 

c 

scfact 

- 

c 

slen 

- 

c 

mm 

- 

c 

velocb 

- 

c 

densb 

- 

c 

dia 

- 

c 

areab 

- 

c 

dzi 

- 

c 

c 

dzr 

- 

c 

dzt 

- 

c 

c 

tstep 

- 

c 

ndp 

— 

c 

tc 

- 

transmitted  wave  stress-time  history 
incident  stress  -  {-  reflected  stress) 
incident  stress  +  ( -  reflected  stress) 
average  specimen  strain-time  history 
average  specimen  strain-rate-time  history 
average  specimen  stress-time  history 
reflected  wave  stress-time  history 
particle  velocity  at  specimen-incident  bar 
interface 

particle  velocity  at  specimen-transmitter  bar 
interface 

interface  stress  difference 
***********  OUTPUT  FILES  *********** 


output  data 

uncorrected  incident  wave  stress-time  history 
uncorrected  reflected  wave  stress-time  history 
uncorrected  transmitted  wave  stress-time 
history 

corrected  incident  wave  stress-time  history 
corrected  reflected  wave  stress-time  history 
corrected  transmitted  wave  stress-time  history 
strain-time  history 
strain-rate-time  history 
stress-time  response 
stress-strain  response 

***********  SYMBOLS  *********** 

gain  for  incident  bar  strain  gauge 

gain  for  transmitter  bar  strain  gauge 

incident  wave  start  time  in  microseconds 

reflected  wave  start  time  in  microseconds 

transmitted  wave  start  time  in  microseconds 

experiment  identification 

output  experiment  identification 

incident  bar  strain  gauge  scale  factor 

transmitter  bar  strain  gauge  scale  factor 

initial  specimen  length 

number  of  Fourier  terms 

pressure  bar  rod  velocity 

density  of  pressure  bar 

diameter  of  pressure  bar 

area  of  pressure  bar 

distance  from  strain  gauge  to 

specimen-incident  bar  interface 

(-)  dzi 

distance  from  strain  gauge  to 
specimen-transmitter  bar  interface 
data  collection  interval 

number  of  data  points  needed  for  computations 
1/tstep 
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# 
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I 

I 


# 
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c  deltat 

c 

c  ibeg 

c  ibegr 

c  ibegt 

c  nlbi 

c  nlbr 

c  nlbt 

c  ibeg i 

c  ibegrr 

c  ibegtt 

c  iend 

c  iendr 

c  iendt 

c  tip 

c  trp 

c  ttp 

c  ni 

c 

c  nr 

c 

c  nt 

c 

c  const 

c 
c 
c 

c  cal  -  computes  scale  factors 

c  basadj  -  zeros  wave  baseline 

c  four  -  performs  dispersion  correction 

c 

c********************************************************** 

C*  MAIN  PROGRAM  * 

c********************************************************** 
c 

program  xx( tty , input=tty,output=tty,cali ,calt, 
lhbi ,  hbt ,  tapell=cal  i ,  tapel2*=calt ,  tapel3=hbi , 
ltapel4=hbt# tape99 ) 
c 

dimension  calt( 4000 ) ,cali ( 4000 ) , bar inc( 25000) , 
ltrsbar( 25000) ,diff (2500) ,sum( 2500) ,xint( 2500) , 
lxrate(2500) ,xstres(2500) ,barref (25000) ,velin(2500) , 
lvelout(2500) , stsdel( 2500) 
c 

1003  format(a) 

1006  format{ fl2.4,lh,fl2.6) 

1007  format('l  ****  hopkinson  bar  data  reduction  code 
loutput  ****•) 

1009  format( 12( lx,f 9 .4) ) 

1010  format('l  ****  processed  hopkinson  bar  data  ****•) 

1011  format(fl0.1,8xff5.1, 3x,f4.1, 3x,f6.1,4x,f 5. 1, 
13x,f4.1,3x,f5.1,2x,f4.1,2x/f5.1) 

character* 20  hopxxi,  hopxxr,  hopxxt,  hopxxif, 


-  2  times  the  wave  travel  time  between  the 
strain  gauge  and  specimen-bar  interface 

-  incident  wave  start  time  in  data 

-  reflected  wave  start  time  in  data 

-  transmitted  wave  start  time  in  data 

-  incident  wave  baseline 

-  reflected  wave  baseline 

-  transmitted  wave  baseline 

-  ibeg  -  nlbi 

-  ibegr  -  nlbr 
-ibegt  -  nlbt 

-  incident  wave  stop  time  in  data 

-  reflected  wave  stop  time  in  data 

-  transmitted  wave  stop  time  in  data 

-  duration  of  incident  wave 

-  duration  of  reflected  wave 

-  duration  of  transmitted  wave 

-  period  of  incident  wave  for  dispersion 
correction 

-  period  of  reflected  wave  for  dispersion 
correction 

-  period  of  transmitted  wave  for  dispersion 
correction 

-  integration  constant 
***********  SUBROUTINES  *********** 


173  - 


lhopxxrf,  hopxxtf , straxx,  strtxx,  strsxx,  sssrxx, 
lexpnam#  hopout 
character*4  test 


operator  supplied  input 


pr int* #' input  experiment  identification  number' 

read  1003# test 

print* 

print*#' gain  for  the  incident  and  transmitter  bar 
lstrain  gauges' 
read*  #gaini #  gaint 
print* 

print*, 'initial  start  time  for  incident  pulse 
1 ( microseconds ) ' 
read*,sttime 
print* 

print* #' initial  specimen  length  (m)' 

read*#slen 

print* 

print* #' input  the  number  of  terms  to  be  used  in' 
print*, 'the  Fourier  analysis' 
read*# mm 

i********************************************************* 

*******  set  file  names  names  for  storage  ******* 

hopxx i= ' hp ' //test// ' i ' 
hopxxr* ' hp ' //test// ' r ' 
hopxxt* ' hp ' //test// ' t ' 
hopxx i f * ' hp ' // tes  t // ' i f ' 
hopxxrf * ' hp ' //test// ' rf  * 
hopxxtf* ' hp ' //test// ' tf ' 
straxx* ' str ' //test// • time ' 
strtxx* ' rat ' //test// ' time ’ 
strsxx* ' sts ' //test// ' time ' 
sssrxx*'prs'//test//' strain' 
expnam* ' hb ' //test// ' x ' 
hopout* ' hp ' //test// ' out ’ 


set  parameter  values 


velocb*4886. 
densb*B090 . 
dia*0. 06033 

areab*( 3.14159*(dia**2))/4.0 

dzi-0.6033 

dzr*-dzi 

dzt*dzr 

tstep*0.5 

ndp*15000 


-  174  - 


wsmm 


s 


tip»240 

trp«tip 

ttp»tip 

delta t«2.*dri/velocb 
srtime»sttime+deltat*le+6 
tstime»srtime 
tc«l/tstep 
c 

c* ********************************************************* 
e  *******  read  in  the  calibration  and  data  files  ******* 
c 

read (11, *) (cali( i) , i«l,4000) 
read (12,*) (calt(i) , i>l,4000) 
read (13,*) (barinc( i) , i*l,ndp) 
read (14,* ) ( trsbar ( i) , i*l,ndp) 
c 

c********************************************************** 

c  ***  calculate  the  scale  factor  from  calibration  data  *** 
c 

call  cal(cali,gaini,scfaci,4000) 
call  cal(calt,gaint,scfact,4000) 
c 

c* ********************************************************* 
c  ****  compute  starting  time  of  waves  in  the  data  **** 
c 

ibeg=tc*sttime+l 

ibegr*tc*srtime+l 

ibegt*tc*tstime+l 

c 

c* ********************************************************* 
c  *****  compute  stop  times  of  waves  in  the  data***** 
c 

iend=tc* ( sttime+tip) +1 
iendr*tc* (srtime+trp)+l 
iendt*tc* ( tstime+ttp)+l 
c 

c* ********************************************************* 
c  *******  length  of  baseline  for  each  wave  ******* 
c 

nlbi-50 

nlbr-20 

nlbt=50 

c 

c* ********************************************************* 
c  ********  compute  baseline  start-times  ******** 

c 

ibeg i«ibeg-nlbi 
ibegrr«ibegr-nlbr 
i beg tt« ibeg t-nlbt 
c 

c*** ******************************************** *********** 
c  *******  tero  wave  baseline  ******** 

c 
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call  basad j( bar inc,scfaci, ibegi, iendr, nlbi, base) 
call  basad j ( trsbar , scfac t , ibeg tt , iendt , nlbt , base ) 
c 

c************************ *************  ******************** 
c  **  creating  a  separate  array  to  store  reflected  pulse  ** 
c 

do  25  i»iend, iendr 
barref ( i)-barinc( i) 

25  continue 
c 

c********************** *********************************** 
c  *******  store  raw  data  ******* 

c 

xi«( ibeg-1) * tstep 
xr*( ibegrr-l)*tstep 
xt*( i beg tt-1)* tstep 
c 

open ( unit*15 , f ile*hopxxi ) 
open ( uni t*16 , f i le=hopxxr ) 
open ( un i t=17 , f i le=hopxxt ) 
c 

do  30  i»ibegiriend 
write(15,1006)xi,barinc( i)/le+6 
xi»xi+tstep 
30  continue 
c 

do  35  i*ibegrr, iendr 
write( 16,1006)xr, barref ( i)/le+6 
xr*xr+tstep 
35  continue 
c 

do  40  i^ibegtt, iendt 
write(17,1006)xt,trsbar( i)/le+6 
xt=xt+tstep 
40  continue 
c 

close (15) 
close (16) 
close (17) 


c  *  perform  dispersion  correction  on  each  wave  separatly  * 
c 

©  nr*( iendr-iend+50) *2 

call  four ( barref , tstep , dia , velocb , nr , dzr , iend , mm) 
c 

ni«( iend-ibegi)*2 

call  four ( bar inc , tstep , dia , velocb , ni , dz i , i beg i , mm) 
c 

^  nt-( iendt-ibegtt) *2 

call  four( trsbar , tstep, dia, velocb ,nt,dzt, ibeg tt, mm) 
c 
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*******  store  dispersion  correction  data  ******* 

42  xi*( ibegi-l)*tstep 
xr»( ibegrr-l)*tstep 
xt»( ibegtt-l)*tstep 
c 

open ( un i t*18 , f i le»hopxxi f ) 
open(unit»19 , f ile»hopxxrf ) 
open ( uni t«20 , f ile*hopxxt f ) 
c 

do  45  i»ibegi,iend 
write (18, 1006 ) xi,bar inc ( i)/le+6 
xi-xi+tstep 

45  continue 
c 

do  47  i«ibegrr , iendr 
write(19,1006)xr,barref ( i)/le+6 
xr*xr+tstep 

47  continue 
c 

do  48  i*ibegtt, iendt 
write ( 20 , 1006 )xt, trsbar( i)/le+6 
xt*xt+tstep 

48  continue 
c 

closet  18) 
close(19) 
close (20) 
c 

c*** *************** ******  *********************  ************* 


c  *****  compute  the  following:  ***** 
c  1)  strain, 

c  2)  strain-rate, 

c  3)  stress, 

c  4)  particle  velocity, 

c  5)  interface  stress  difference 

c 


j«ibegr 

ib*l 

do  50  i*ibeg,iend 
diff ( ib)=barinc( i)-barref ( j) 
sum(  ib)**barinc(  i)+barref  (  j) 
j-j+1 
ib»ib+l 
50  continue 
c 

const*l/( slen*densb*velocb) 
c 

ipoint*iend-ibeg 

xint(l)«0.0 

xrate(l)«0.0 

tstepi«tstep*1.0e-6 

ic-2 
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j-ibegt+1 

do  100  i>ictipoint 
n«i-l 

k-j-1 

xratet i)»(dif f (n)-trsbar (k) )*const 

xint( i)«xint(n)+( tstepi* 0.5* (xratet i ) +xrate ( n ) ) ) 

j-j+1 

100  continue 
c 

j«ibegt 

do  120  i«l,ipoint 
xstres( i)*(sum( i)+trsbar( j) )/2.0 
j«j+l 

120  continue 
c 

do  127  i«l,ipoint 
n*i-l 

velin( i)»( (dif f ( i) )/(densb*velocb) ) 

velout ( i)«( ( trsbar ( ibegt+n) )/(densb*velocb) ) 

127  continue 
c 

do  128  isl,ipoint 
n*i-l 

stsdelt i)*( sumt i ) -trsbar ( ibegt+n) )/le+6 

128  continue 
c 

c* ******************************************************** 
*******  write  output  files  ******* 

open ( uni t*21,f ile»straxx) 
open(unit*22 ,f ile=strtxx) 
open(unit»23/f ile*strsxx) 
open(unit=24,f ile*sssrxx) 
t*0.0 

do  130  i«l,ipoint 
xstres ( i )=xstres ( i)/le+6 
xint( i)*xint( i) *100. 
write(21,1006)t,xint< i) 
wr i te ( 22, 1006 )t, xratet i) 
write (2 3,10 06)t, xstres( i) 
write(24,1006)xint( i) , xstres ( i) 
t»t+tstep 
130  continue 
close(21) 
closet  22) 
close(23) 
closet  24) 


c 

c 

c 

c 


*******  output 

xi-( ibeg-l)*tstep 
xr*( ibegr-1) *tstep 
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c 


c 


c 


xt»{ ibegt-1 ) *tstep 


open ( un i t*2  7 , f i le=hopout ) 
write(27,1007) 


writet 27, 
write(27, 
writet 27 , 
writet 27/ 
writet  27 , 
write(27, 
writet 27, 

1 - 

writet 27, 
write(27, 

1 - 

writet 27, 
writet 27 , 

1 -  ',g 

writet 27, 
writet 27, 

1 -  '  ,xi 

writet27, 
writet  27, 

1 -  '  ,xr 

writet  27, 
writet  27 , 
1-  * » xt 

writet  27 , 
writ-t  27 , 

1 - 

writet  27, 
writet27, 

1 —  1  ,  mm 

writet27, 
writet  27 , 

1 - 

writet  27 , 
writet 27, 

1 - 

writet  27, 
writet27, 

1 - 


***  input  parameters  ***• 


\  t 
\  I 
\  I 
\  I 
\  I 
\  I 

) 'experimental  identification  number 
'  , expnam 
)  1  1 

)  'gain  for  incident  bar  strain  gauge 
' ,gaini 
) '  ' 

)'gaint  for  transmitter  bar  strain  gauge 
int 

\  I  I 

) 'incident  wave  starttime  t microseconds) 


wave  starttime  t microseconds) 


)'  ' 

) 'reflected 
)'  ' 

) 'transmitted  wave  starttime  t microseconds) 
)'  ’ 

) 'initial  specimen  length  (m) 

-  ',slen 

)'  ' 

) 'number  of  terms  used  in  fourier  analysis 


)'  ' 

) 'dispersion  bar  correction  length(m) 
'  ,dzi 
)'  ' 

) 'incident  bar  calibration  factor 
' , scf aci 
)’  ' 

) 'transmitter  bar  calibration  factor 
', scf act 


write(27,1010) 
write( 27 ,* ) '  ' 
writet 27,*)'  ' 
writet 27,*)'  ' 

writet27,*)'  time  stress  strain  strainrate 

linfll  veil  inf 12  vel2  stsdif* 
writet 27 ,*)' tmicroseconds )  tmpa)  t%)  tl/sec) 

l(mpa)  (m/s)  tmpa)  (m/s)  tmpa)' 
write(27,*)'  ' 
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ti*0 

do  137  i>l,ipoint 
n-i-1 

write ( 27 , 1011 ) ti ,xstres( i) ,xint( i) ,xrate( i) , 
lsum( i)/le+6, velin( i) , trsbar( ibegt+n)/le+6, 
lvelout(  i)  fStsdeK  i) 
ti»ti+tstep 
137  continue 

close (27) 

4  end 


END  MAIN  PROGRAM  * 

c* ********************************************************* 
C 

c 

c********************************************************** 
C*  SUBROUTINES  * 

c***«* ***************************************** ************ 
C 

c 

c  This  subroutine  is  used  to  compute  scale  factors 

c  that  are  used  to  transform  the  digitized  numbers 

c  to  stress-time  data 
c 
c 

subroutine  cal ( blk ,  gain , self ac , ndp ) 

dimension  blk(4000) 

calmax*0.0 

calzer»0.0 

ncalmx=0 

ncalze=0 

ibeg«l 

if (blk(l) .gt .164)  goto  130 
100  do  105  i*ibeg+l,ndp 

if (blk( i) .gt.164)  then 
iend^i 
goto  110 
endif 

105  continue 
iend«ndp 

110  if ( (iend-ibeg) .gt.8)  then 
do  115  i»ibeg+4, iend-4 
calzer«calzer+blk( i) 

115  continue 

ncalze*ncalze+iend-ibeg-7 

endif 

ibeg>iend 

if ( ibeg.ge.ndp)  goto  190 
130  do  135  i«ibeg+lf ndp 


u  o 


4ft 


ft 


i£(bl)c( i)  .It. 164)  then 
iend«i 
goto  140 
endif 

135  continue 
iend^ndp 

140  if( ( iend-ibeg ) .gt . 8 )  then 
do  145  i«ibeg+4, iend-4 
calmax=calmax+blk( i) 

145  continue 

ncalmxancalmx+iend-ibeg-7 

endif 

ibeg»iend 

if ( ibeg.ge.ndp)  goto  190 
goto  100 

190  calmax«calmax/ncalmx 
calzer*calzer/ncalze 
self ac=l . 2611e+9/gain/( calmax-calzer ) 

return 

end 


# 


I 

i 


e 


» 


c 

c 

c 

c 

c 


This  subroutine  perform  the  dispersion  correction 


subroutine  four( fet, tstep,dia,velocb,nn,dz, 
libegin,mm) 


c 

c 

c 

c 

fct( i) 

c 

nn 

c 

a 

c 

b 

c 

t 

c 

c 

c 

c 

equa  ( i ) 

c 

V 

c 

wo 

c 

ao 

c 

d 

c 

c 

CO 

c 

c 

c 

c 

coef  f 

*********  SYMBOLS  ********** 

■  array  holding  wave  data 

■  period  of  the  wave  being  analyzed 

*  array  to  store  An  Fourier  coefficient  values 

*  array  to  store  Bn  Fourier  coefficient  values 

*  array  to  store  phase  angle  values 

■  array  that  holds  the  phase  velocity  at  each 
frequency 

■  dummy  array  for  computational  purposes 

■  bar  radius/wavelength 

■  wave  frequency 

*  value  of  Ao  Fourier  coefficient 

■  array  to  store  the  value  of  Cn  after 
transformation  to  a  cosine  series 

*  longitudinal  wave  velocity  of  bar 

*******  SUBROUTINES  ******* 

-  computes  the  Fourier  coefficients  Ao,  An, 
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I 


* 


C 


c  and  Bn 

c  simp  -  performs  Simpson*  rule  integration 

c 

c 

dimension  fct(25000) ,a(100) ,b(100) ,c(100) , 
ld(100) ,t(100) 
c 

c  function  statement;  phase  velocity  verses  r/lambda 
c 

cp(v)-.5764+( . 4236/(22 .*v** 4+12 .8*v**3-2 .77*v**2+ 
1.92*v**l. 5+1.0)) 
c 

f req-1 .0/1 tstep*1.0e-6) 
c 

do  20  j=l,nn 
n*j-l 

fct( j)«fct( ibegin+n) 

20  continue 
pi-3.14159 
wo-2*pi/nn 
c 

call  coeff (mm,nn,wo,fct ,ao,a,b) 
c 

do  60  j*l ,mm 

d( j)-sqrt(a( j) **2+b( j)**2) 
t( j)-acos(a( j)/d( j ) ) 
if (b( j)-0.0)  55,60,60 
55  t(  j)—  1.0*t(  j) 

60  continue 
c 

69  w-(dia/2.0)*wo*freq/velocb 
v»0.0 
i-0 

70  i-i+1 
x=i*w 

80  v-v+0.0001 

xp-2.0*pi*cp< v)*v 
83  if (xp-x)  80,85,85 
85  cCi)-cp(v) 

88  if (i-mm)70,90,90 

90  do  100  i-l,mm 

c( i)-c( i)*velocb 
100  continue 

do  120  i-l,mm 

ak»  i  *wo/velocb*  ( velocb/c  ( i  >  -1 . 0 )  *‘f  req 
t( i)-t( i)+ak*dz 
120  continue 

perform  wave  reconstruction 
3  do  150  i-l,nn 
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e 


fct{ i)«ao/2.0 
do  150  j“l ,mm 

fct(i)-fct{i)+d( j)*cos( ( j*wo*i)-t( j) ) 

ISO  continue 
c 

do  160  j>l,nn 
n*  j-1 

fct( ibegin+n)«fct( j) 

160  continue 
c 

return 

end 

c 

c 

c*** **************************************************** 

c 

c 

subroutine  coeff  tnun,nn,wo, fct,ao,a,b) 
c 

c  this  routine  will  compute  the  Fourier  coefficients 
c  for  each  term  in  the  series 
c 

dimension  equal 6000) , a(100) ,b(100) ,fct( 25000) 
c 

c  computing  the  integration  limit;  t/2 
c 

n=nn/2 

do  10  i«l,nn 
equal i)»fctt i) 

10  continue 
c 

c  Simpson's  rule  integration  will  be  performed  to 

c  evaluate  each  of  the  coefficients;  Ao,An,Bn 

c 

c  computing  Ao 

c 

call  simp(equa,nn, result) 
c 

ao»result/n 

c 

c  compute  An 

c 

do  30  j*l ,nun 
do  20  i*i,nn 

equal i)»fct( i)*co8l i*j*wo) 

20  continue 

call  simplequa,nn, result) 
at j)«result/n 
30  continue 

compute  Bn 


do  50 
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do  40  1*1, nn 

equal i)»fct( i)*sin( i* j*wo) 

40  continue 

call  s imp (equa,nn, result) 
b( j)*result/n 
50  continue 
c 

return 

end 

c 

c 

c******************************************************** 

c 

c 

subroutine  s imp (equa,nn, result) 
dimension  £(6000) ,equa(6000) 
c 

do  2  i»l,nn 
f (i)-equa(i) 

2  continue 
c 

n«nn/2 

c 

c  this  subroutine  performs  Simpson's  rule  integration 
c 

c  parameters 

c 

c  £:  integration  points 

c  n:  number  of  data  points 

c  h:  spacing  between  data  points 

cresult:  estimate  of  integral  returned  to  caller 

c 

c  check  to  see  if  the  number  of  panels  is  even 
c 

c  n-1:  number  of  panels 

c 

npanel-n-1 

nhalf-npanel/2 

nbegin*l 

result-0.0 

h-1 

if ( (npanel-2*nhalf ) .eq.0)  go  to  5 
c 

c  if  the  number  of  panels  is  odd  use  3/8  rule  for  the 
c  first  four  points  then  continue  using  the  1/3  rule 

c  if  the  numbers  of  panels  is  even  then  the  1/3  rule 

c  will  be  used  exclusively 
c 

result-3. 0*h/8.0*(f <l)+3.0*f (2)+3.0*f (3)+f (4) ) 
nbegin-4 

apply  1/3  rule,  add  in  the  first, second,  and  last 
values 


nooon  o  non 


c 


5  result-resul t+h/3 . 0* ( £ ( nbeg in )  + 

14. 0*f (nbegin+1 )+f ( n) ) 
nbegin«nbegin+2 
if (nbegin.eq.n)  go  to  15 

the  pattern  after  nbegin+2  is  repetitive 

nend«n-2 

do  10  i>nbegin,nend,2 
result»result+h/3.0*(2.0*f ( i)+4.0*f ( i+1) ) 
10  continue 

15  return 
end 


********************************************************* 


i  • 


» 


C 


subroutine  basad j ( blk , self ac , nbeg , nend , n , base ) 
dimension  blk( 25000) 
c 

50  k»nbeg+0.8*float(n) 
base»0.0 
count»0.0 
do  100  i*nbeg,k 
count*count+l 
base*base+blk( i) 

100  continue 
c 

base*base/count 

c 

c  adjust  data  for  baseline  and  apply  scalefactor 
c 

150  do  200  i*nbeg,nend 

blk( i)-(blk( i) -base) *sclf ac 
200  continue 
c 

return 

end 

c 

c 

c*** ****** ****************************************** 
c*************************************************** 
C*  END 

c**************************** *********************** 


* 


) 
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G.3  Computer  Output  for  Experiment  134 

The  following  abbreviations  have  been  used  in  the 
computer  output  (see  also  figure  2.5): 
infll  :  stress  at  interface  1, 
veil  :  particle  velocity  at  interface  1 
infl2  :  stress  at  interface  2 
vel2  :  particle  velocity  ar  interface  2 
stsdif  :  stress  difference  between  interface  1  and 


interface  2 


XXXX  hopkinson  bar  data  raduetion  coda  output  KXXK 


XXX  input  paramatars  XXX 

axparimantal  idantifieation  number  - - -  hbl34ax 

gain  for  incident  bar  strain  gauge  - - —  2.5 

gaint  for  transmitter  bar  strain  gauge  -  5. 

incident  wave  starttime  (microseconds)  -  2519. 

reflected  wave  starttime  (microseconds)  -  2765.5 

transmitted  wave  starttime  (microseconds)  -  2765.5 

initial  specimen  length  (m)  - —  1.26e-2 

number  of  terms  used  in  fourier  analysis  —  20 

dispersion  bar  corraetion  length(m)  - -  0.6033 

incident  bar  calibration  factor  -  7630275.099424 

transmitter  bar  calibration  factor  - - —  3816424.365323 
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APPENDIX  H 

SPECIMEN  MOISTURE/DENSITY  RELATIONSHIPS 
TO  THE  HARVARD  MINIATURE 
COMPACTION  CURVE 


jr'J? 


Figure  H.3.  Moisture/density  relationship  to  the  Harvard 
miniature  compaction  curve  for  specimens 
compacted  wet  of  optimum. 


APPENDIX  I 


PLOTTED  EXPERIMENTAL  RESULTS 
WITH  SUMMARY  DATA 


The  following  nomenclature  has  been  used  in  the 
presentation  of  the  summary  data: 

oa  *  magnitude  of  applied  stress  wave, 

C max  m  maximum  average  strain-rate, 

Cavg  *  average  strain-rate  over  the  duration  of  the 
experiment, 

1  «  specimen  length  after  seating  has  been  applied, 
Cs«  seating  strain, 

w  *  initial  specimen  moisture  content, 

W|  ■  percentage  if  initial  moisture  lost  during  the 
experiment, 

e  ■  void  ratio  based  on  specimen  length  1, 
ng  »  initial  gas  porosity, 

Vt  ■  specimen  total  volume, 

7t  »  wet  density. 


S  ■  degree  of  saturation 


Figure  1.1.  Stress-strain  response  for  experiment  31 
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Figure  1.3.  Stress-strain  response  for  experiment  38 


Figure  1.4.  Stress-strain  response  for  experiment  39 


Stress-strain  response  for  experiment  57 
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Stress-strain  response  for  experiment  61 


270  MPa 


C»M 


Figure  1.28.  Stress-strain  response  for  experiment  118 


375  MP 


Stress-strain  response  for  experiment  138 


523  MPa 


237  MPa 


522  MPa 


507  MPa 
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Due  to  the  uncertainty  in  the  evaluation  of  the 
strain  contributions  from  factors  discussed  in  section 
6.4,  a  gross  strain  adjustment,  to  account  for  the  excess 
between  strain  at  peak  stress,  and  the  initial  gas 
porosity  has  been  computed.  This  was  accomplished  by 
plotting  the  difference  between  the  strain  at  peak  stress 
and  the  initial  gas  porosity  against  the  average  force 
felt  by  the  specimen.  The  average  force  is  computed  as  the 
peak  stress  multiplied  by  the  area  of  the  specimen.  A 
linear  regression  line  was  then  fit  to  the  data  for 
specimens  prepared  near  the  optimum  compaction  conditions. 
The  linear  regression  line  was  computed  with  the  dependent 
variable  (Y)  taken  as  the  difference  between  the  initial 
gas  porosity  and  the  strain  at  peak  stress  and  the 
independent  variable  (X)  taken  as  the  maximum  force  felt 
by  the  specimen.  The  equation  of  the  linear  regression 
line  is; 

Y  «  0.9BX  +  0.56  .  (J.l) 

The  resulting  plot  is  shown  in  figure  J.l.  The  data  used 
to  prepare  figure  J.l  are  presented  in  table  J.l. 

Using  the  regression  line  and  the  average  force 
sustained  by  the  specimen,  a  strain  correction  was 
computed  for  each  specimen.  The  strain  correction  is  added 
to  the  initial  gas  porosity  and  the  sum  compared  to  the 
strain  at  peak  stress.  The  balance  can  then  be  used  to 


Table  J.l 


Strain  Correction  Data 


Experiment  No. 

Strain  -  Gas  Porosity 
(«) 

Force 

(dynes  10E10) 

112 

11.9 

7.5 

113 

9.7 

7.2 

114 

10.4 

7.2 

115 

7.5 

7.7 

116 

6.2 

4.5 

117 

6.0 

6.1 

118 

3.8 

4.7 

119 

4.2 

4.5 

131 

3.1 

5.6 

132 

3.8 

6.0 

133 

4.5 

5.1 

134 

6.4 

6.1 

138 

6.9 

5.3 

145 

3.9 

3.0 

146 

7.0 

9.0 

147 

6.3 

4.3 

162 

16.2 

9.4 

163 

10.4 

12.2 

164 

4.9 

8.9 

165 

2.2 

2.2 

a)  Strain  is  taken  as  the  strain  at  peak  stress. 
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gauge  the  success  of  the  correction.  The  results  of  the 
correction  procedure  are  tabulated  in  table  J.2. 


This  procedure  requires  none  of  the  assumptions  made 
in  the  computations  of  strain  contribution  from  moisture 
loss,  pore  water  compression,  or  radial  expansion.  It  is 
based  solely  on  the  peak  stress  sustained  by  the  specimen. 
In  addition,  the  strain  contribution  of  soil  loss  is  also 
included. 


* 


Table  J.2 


Strain  Correction  Results 


Experi—nt 

Cm  Porosity 

Strain 

Bub 

Strain  at 

Balance 

MO. 

<») 

Correction 

Paak  Stress 

(•) 

(«) 

(«) 

(1) 

• 

112 

7.91 

7.87 

15.78 

19.79 

4.01 

113 

7.73 

7.67 

15.40 

17.40 

2.00 

114 

9.a 

7.61 

16.82 

19.60 

2.78 

115 

10.55 

8.07 

18.62 

18.03 

-0.59 

116 

8.07 

4.98 

13.05 

14.24 

1.19 

117 

8.50 

6.54 

15.04 

14.50 

-0.54 

118 

10.25 

5.18 

15.43 

14.06 

-1.37 

• 

119 

9.91 

4.93 

14.84 

14.07 

-0.77 

131 

9.57 

6.10 

15.67 

12.68 

-2.99 

132 

6.88 

6.43 

13.31 

10.67 

-2.64 

133 

8.02 

5.52 

13.54 

12.52 

-1.02 

134 

S.80 

6.57 

12.37 

12.22 

-0.15 

138 

4.86 

5.75 

10.61 

11.80 

1.19 

145 

6.23 

3.51 

9.74 

10.12 

0.38 

146 

6.26 

9.41 

15.67 

13.30 

-2.37 

• 

147 

6.64 

4.78 

11.42 

12.99 

1.57 

162 

4.44 

9.74 

14.18 

20.60 

6.42 

163 

4.84 

12.52 

17.36 

15.27 

-2.09 

164 

7.49 

9.28 

16.77 

12.43 

-4.34 

165 

7.05 

2.70 

9.75 

9.27 

-0.48 

• 

# 

e 

i 

• 
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